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Abstract. It is well known that subspaces of the Hardy space over the unit 
disk which are invariant under the backward shift occur as the image of an ob- 
servability operator associated with a discrete-time linear system with stable 
state-dynamics, as well as the functional-model space for a Hilbert space con- 
traction operator. We discuss two multivariable extensions of this structure, 
where the classical Hardy space is replaced by (1) the Fock space of formal 
power series in a collection of d noncommuting indeterminates with norm- 
square-summable vector coefficients, and (2) the reproducing kernel Hilbert 
space (often now called the Arveson space) over the unit ball in C d with repro- 
ducing kernel fc(A,C) = 1/(1- (A, 0) (\c G C d with ||A||, HCII < 1). In the first 
case, the associated linear system is of noncommutative Fornasini-Marchesini 
type with evolution along a free semigroup with d generators, while in the sec- 
ond case the linear system is a standard (commutative) Fornasini-Marchesini- 
type system with evolution along the integer lattice Z d . An abelianization 
map (or symmetrization of the Fock space) links the first case with the second. 
The second case has special features depending on whether the operator-tuple 
defining the state dynamics is commutative or not. The paper focuses on 
multidimensional state-output linear systems and the associated observability 
operators; followup papers 9 10'j use the results here to extend the analysis 
to represent observability-operator ranges as reproducing kernel Hilbert spaces 
with reproducing kernels constructed from the transfer function of a conserva- 
tive multidimensional (noncommutative or commutative) input-state-output 
linear system. 
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1. Introduction 

For U and y any pair of Hilbert spaces, we use the notation C(U, y) to denote 
the space of bounded, linear operators from U to y. For X a single Hilbert space, 
we shorten the notation C(X, X) to C(X). Let X, U and y be Hilbert spaces, 
let A e C(X), B e C(U,X), C e C(X,y) and D e £(U,y) be bounded linear 
operators, and let us consider the associated discrete-time linear time-invariant 
system 

x(n + 1) = Ax(n) + Bu(n) (11) 
y{ri) = Cx(n) + Du(n) 

with x(n) taking values in the state space X, u(n) taking values in the input-space 
IA and y(n) taking values in the output-space y. If we let the system evolve on the 
nonnegative integers n S Z+, then the whole trajectory {u(n), x(n), y(n)} n <£i + is 
determined from the input signal {u(n)} n £% + and the initial state x(0) according 
to the formulas 

n-1 

x(n) = A n x{0) + A^^Buik), 

k=0 
n-1 

y(n) = CA n x(0) + ^ CA n - l - k Bu{k) + Du{n) 

k=0 

n-1 

= [O c ,Ax(Q)] n + CA"- 1 - k Bu(k) + Du(n) (1.2) 

k=0 

where Oc,a denotes the so-called observability operator 

Oc, A :x^{CA n x} n&+ . 
If we introduce the Z-transform 

oo 

U(n)} n&+ - /(A) = /(") A "' (1-3) 

the Z-transformed version of the system-trajectory formulas (|1.2|) become 

x{\) = {I- XAy 1 x(0) + A(7 - \A)~ 1 Bu{X), 

y(X) = C(I - Ai4) -1 a;(0) + [D + XC(I - AA)- x B]u(A) 

= a c , A a:(0)+r E (z)u(A) (1.4) 

where 

oo 

O c ,A- x^^{CA n x)X n =C(I-XA)- 1 x (1.5) 
n=0 

is the Z-transformed version of the observability operator and where 

T S (A) = D + XC(I - XA^B 

is the transfer function of the system S given by In particular, if the input 

signal {ii(n)}„ e z + is taken to be zero, the resulting output {y(n)} n ^z + is given 
by y = Oc,ax(0). In case Oc,a is bounded as an operator from X into i y :— 
1 2 <g> y (here f 2 is the space of square-summable complex sequences indexed by 
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the nonnegative integers Z+ , we say that the pair (C, A) is output- stable. It is 
convenient to represent Oc,A in the output-stable case in the matrix form 

O c ,a = col„ 6Z+ [CA n ] : X^l%, 

Since the Z -transform maps £y unitarily onto Hy := H 2 (H)) ® 3^, where i? 2 , 
the image of £ 2 under the ^-transform, is the space of analytic functions on the 
unit disk with modulus-square-summable sequence of Taylor coefficients: 

oo oo 

h 2 = {/(A) = Yl /« A " : E i/«i 2 < °°}< 

n=0 n=0 

the output stability of (C, A) is equivalent to the Z-transformed version of the 
observability operator (|1.5|l being bounded as an operator from X into Hy, It is 

readily seen that Ocax — Oc,ax. 

If (C, A) is output-stable, then the observability gramian 

Qc,a ■= (Oc,a)*Oc,a = {Oc,a)*Oc,a 
is bounded on X and can be represented via the series 

g c ,A = Y, A * nc * CAn ( L6 ) 

71=0 

converging in the strong operator topology. The following result gives a summary 
of well-known connections between output stability, observability gramians and 
solutions of associated Stein equations and inequalities. 

Theorem 1.1. Let (C, A) be a pair of operators with C: X — > y and A: X — > X . 
Then: 

(1) The pair (C, A) is output-stable if and only if the Stein inequality 

H - A* HA > C*C (1.7) 

has a positive semidefinite solution H 6 C(X). 

(2) If (C, A) is output- stable, then the observability gramian Qc.A satisfies the 
Stein equality 

H - A* HA = C*C (1.8) 
and is the minimal positive semidefinite solution of the Stein inequality 

O- 

(3) There is a unique positive semidefinite solution of the Stein equality (11.8(1 
if A is strongly stable, i.e., powers A n of A tend to zero in the strong 
operator topology of C{X). If A is a contraction operator, then the positive 
semidefinite solution of the Stein equation (|1.8|l is unique if and only if A 
is strongly stable. 

A pair (C, A) is called observable if the operator Oc,a (equivalently, Oc,a, Sc,a) 
is injective. This property means that a state space vector x G X is uniquely 
recovered from the output string {y(n)}^L generated by running the system 
with the zero input string and the initial condition x(0) = x. A pair (C, A) is called 
exactly observable if Oc.a (equivalently, Gc.a) is bounded and bounded from below. 

Associated with an output-stable pair (C, A) is the range of the observability 
operator 

Ran 6 C ,A = {C(I - zA)~ x x : x G X}. 
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The following theorem summarizes the connection between such ranges and back- 
ward-shift-invariant subspaces of Hy. 

Theorem 1.2. Suppose that (C, A) is an output-stable pair. Then: 

(1) The linear manifold KanOc.A is invariant under the backward shift opera- 
tor 

S':/«.M. (1.9, 

(2) Let H > be a solution of the Stein inequality 11. 7f) and let X' be the 
completion of X with inner product \\[x\\\ 2 X i = {Hx,x)x (where [x] denotes 
the equivalence class modulo Keri/ generated by x). Then A and C extend 
to define bounded operators A' : X' — * X' and C ' : X' — > y and the observ- 
ability operator Oc,a extends to define a contraction operator Oc' t A' from 
X' into Hy. Moreover, Oc'.A' ■ X 1 — > Hy is an isometry if and only if H 
satisfies the Stein equation p. 8(1 and A' is strongly stable, i.e., 

(HA n x, A n x) for all x € X. 

(3) // the linear manifold Ai := Ran Oq .a is given the lifted norm 

0c,a x \\m = ^ v „ M „ ( H y>y)x, 

then 

(a) M can be completed to M! = RanO^*/ # with contractive inclusion in 
H y : 

ll/llffj < UWm' for allfeM'. 
Furthermore, Ai' is isometrically equal to the reproducing kernel Hil- 
bert space with reproducing kernel Kc,a,h given by 

K c ,aM\ = C(I - XA^Hil - (A*)- 1 ^. (1.10) 

(b) The following difference- quotient inequality is satisfied 

\\S*f\\ 2 M <\\f\\ 2 M-\\f(0)\\y for allfeM (1.11) 
and moreover, if the Stein equality (|1.8|) holds, then Ijl.llfl holds with 
equality. 

(4) Conversely, if Ai is a Hilbert space contractively included in Hy which is 
invariant under S* and for which the difference- quotient inequality (|1.11|) 
holds, then there is a contractive pair (C, A) (i.e., I|1.7I) holds with H = 1%) 
such that Ai = TC(Kc.aj) = RanOc^ isometrically. In case HI- 11( 1 holds 
with equality, then (C, A) can be taken to be isometric. 

Results of the type in Theorem 1 1 . 1 1 are the basis for the Lyapunov- function ap- 
proach to stability analysis in system theory; there are far-reaching generalizations 
to nonlinear and time-varying settings which are far afield from our main interests 
here. The goal of characterizing subspaces of H 2 of the form TC(Kc,a-,h) (especially 
in a finite-dimensional context) was a key feature in the approach to Nevanlinna- 
Pick interpolation developed by H. Dym (see |29p. 

In this paper we present the analogues of Theorems 11.11 and 11.21 for the two 
related multivariable settings: (1) the case where the Hardy space H 2 on the unit 
disk is replaced by the Fock space Hy(Td), and (2) the case where H 2 is replaced 
by the vector- valued Arveson reproducing kernel Hilbert space Hy(kd). 
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To define the Fock space, we let T d denote the free semigroup on the set {1, . . . , d} 
of the first d natural numbers and then let Hy (Td) consist of the space of all formal 
power series $^ 0& p f yZV m ^ noncommuting indeterminates z = (z%, . . . , Zd) with 
coefficients f v in a coefficient Hilbert space y which are square-summable in norm: 

^v£? d WfAy < °°. Here we write z " = z i N z iN-i ■ ■ - z ii ir v = Ww-i ...ii e T d - 
The shift operator S: f(\) i— > A/(A) acting on the Hardy space H 2 is replaced by 
the noncommuting d-tuple S = (Si, . . . , S d ) on Hy(T d ) given by 

Sj:f(z)~f(z)zj for j = l,...,d. (1.12) 

The system is replaced by a noncommutative multidimensional input-state- 

output system of the form 

{x(lv) = Axx(v) + Bxu(v) 
'■ '■ (1.13) 
x(dv) = A d x(v) + B d u(v) 
y(v) = Cx(v) +Du(v). 

Here the system evolves along the free semigroup and, for each v G Td, the 
state vector x(v), input signal u(v) and output signal y(y) take values in the state 
space X , input space U and output space y, and the system matrix U has the form 











'Ax 


Bi 




'X' 
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B 










X 






C 


D 




A d 


B d 




u 




X 










C 


D _ 




y 



(1.14) 



Such systems were introduced in [THj and with further elaboration in JI] and ; 
following jll| we call this type of system a noncommutative Fornasini-Marchesini 
linear system. The observability operator associated with an output map C : X — > 3^ 
and a d-tuple A = (A\, . . . ,A d ) of not necessarily commuting operators on a Hilbert 
space X , expressed in "frequency-domain" coordinates, takes the form 

O c ,A :ihC(/- ziAx z d A d )~ l x. 

For the particular case where A is a row contraction and 

C = {I-A* X A X A d A d yl 2 

with y taken to be equal to the closure of the range of C, this operator appears 
already in work of Popescu 0Hj under the term "Poisson kernel" and as the adjoint 
of the key operator L used in many constructions in the paper of Arveson 
Reproducing kernel Hilbert spaces consisting of formal power series were developed 
in a systematic way in 15 . Such spaces already appear (although not quite in 
our notation) in the Sz.-Nagy-Foia§ model theory for row contractions developed 
by Popescu (see ^]021E31)- We shall see that Theorems 11.11 and 11.21 extend in a 
natural way to this setting, where the observability gramian i|1.6f) in the statement 
of Theorem 1 1.1 1 is replaced with the multivariable observability gramian 

Sc.a = 6c,a<5c,a = ]T A V *C*CA V (1.15) 

(here we set A" — Ai N ■ ■ ■ A^ if v = iff ■ ■ -ix € T d ), where the backward shift S* 
(|1.9(l in the statement of Theorem 1 1.21 is replaced by the d-tuple S* = (S^, . . . ,S d ) 
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of adjoints of the shift operators Sj in (|1.12(l . and where the positive kernel (|1.1Q(I 
becomes the kernel 

K c ,am(z, w) = C(I - z x A x z d A d )- l H{I - w\A\ waA^C (1.16) 

in two sets z = [z%, . . . , z d ) and w = (lOi, . . . , Wd) of noncommuting indeterminates 
(see Theorems E21 and EH below). 

In the second Arveson-space setting, the Hardy space H 2 over the unit disk is 
replaced by the so-called Arveson space, the reproducing kernel Hilbert space H(kd) 
over the unit ball M d in complex d-dimensional space C d based on the reproducing 
kernel function 

= — ^-=7— for A, c e v d 

and the classical Hardy-space shift /(A) i— ► A • /(A) is replaced by the d-tuple of 
Arveson shift operators Ma = (-^Ai > • • • j M\ d ) where 

M A , : /(A) -> \ j f(\) for / 6 H(k d ) (1.17) 

(see [2E1E])- In this case the underlying system evolves along the integer lattice 
U\_ = (n = (m, • • • , rid) : rij e Z+} has the form of what we call a (commutative) 
Fornasini-Marchesini system 

x(n) = Aix(n-ei)H h A d x(n - e d ) 

+Biu{ti - ei) H h B d u(n - e d ) (1.18) 

y(n) = Cx(n) + Du(n). 

Here and in what follows, ej denotes the element in Zl having the j-th partial 
index equal to one and all other partial indices equal to zero: 

ej = (0, 0,1,0, ...,0) e 1%. (1.19) 

Thus the system matrix U has the same form (I1.14f) as for the noncommutative 
setting but the domain for all the signals and the system evolution is the inte- 
ger lattice rather than the free semigroup Ta. and the associated "frequency- 
domain" objects are functions or formal power series in the commuting variables 
A = (Ai, . . . , Xd) rather than in the noncommuting indeterminates z — (zi, . . . , z d ). 
In Section |3| we show how the Arveson space Tt(kd) ® y and this Fornasini- 
Marchesini linear system can be derived as an abelianization (sometimes also called 
symmetrization) of the noncommutative Fock space H^Td) and of the noncommu- 
tative Fornasini-Marchesini linear system, respectively; while it is well known that 
the Arveson space is a symmetrization of the Fock space and that the multiplier 
algebra on the Arveson space is the image under a completely positive map acting 
on the noncommutative multiplier algebra on the Fock space (see [7] ES| and 
|491 150| for a recent, more general systematic framework), our extension of these 
ideas to the underlying system theory appears to be new. The observability oper- 
ator, as in the noncommutative setting, is associated with a so-called output pair 
(C, A) but now has the form 

Oc,a :x~C(I- XiAt XdAd)-^ 

where the variables Ai, . . . , Xd commute and the abelianized observability gramian 
Qq a = (Oq a )*Oq A has an infinite-series representation more complicated than 
the second expression in (|1.15l) (see equation (|3.10|) below). In case the operator 
rf-tuplc A = (A\, . . . , Ad) is commutative (so AiAj — AjAi for all 1 < i,j < d), 
&c,a = Gc,A (see Proposition !^. HI below), and Theorem 1 1.21 has a natural analogue 
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for this setting, with the abelianized multivariable observability gramian Qq a = 
Gc.a playing the role of the observability gramian in Theorem ll.il with the operator 
d-tuplc = (M£ , . . . , M% d ) in place of the backward shift S* in Theorem 

Ol and with kernel (|0)|) now taken to be the multivariable positive kernel 



Kc,A;h(*> = C(I - XxAx \ d A d )- l H{I - CiA* QA^C* (1.20) 



(see Theorems l3.14ll3.15l and r3.16l below'). In the general case where the operator d- 
tuple A = (Ax, ■ ■ ■ , A d ) is not assumed to be commutative, there is no characteriza- 
tion of the abelianized observability gramian as a minimal solution of a generalized 
Stein equation analogous to the classical case given in Theorem ll.il but there still 
is a somewhat more implicit analogue of Theorem 11.21 where the backward shift 
S* (|1.9|) in Theorem ll.2l is replaced by a solution of the so-called Gleason problem 
(see Theorems 13 .201 and 13 . 2 II below) . The Gleason problem originates in the work 
of Henkin and Gleason (see EBJ) and has been studied in the context of the 
Arveson space (with various formulas for the solution) in jjjj with an application 
to realization questions in |SJ. Our analogue of Theorem 1 1.21 for the Arveson space 
for the case of commutative <i-tuple A has already been given in |19| (with a more 
general power-series setting worked out in |20j) for the finite-dimensional case. 

We also give various numerical examples (constructed with the aid of the soft- 
ware program MATHEMATICA) to illustrate how Oc,a an d Oq a can have divergent 
properties when A is not commutative (see Examples 13.41 l3~9*l and |3~TD below). 

Backward-shift-invariant subspaces for the classical setting have been used for 
some time in the operator-theory literature as the model space for a more general 
(abstractly defined) Hilbert-space contraction operator (see [21001); connections 
of this work with linear system theory were only realized later (see e.g. |34l 1351 ). 
Our results develop the structure of such model spaces for the case of operator- 
tuples and therefore are of interest from the point of view of multivariable operator 
theory. We find it satisfying that these model spaces in turn tie in with the theory 
of multidimensional linear systems in much the same way (but with some surprises) 
as in the classical case. 

As applications of the ideas, we obtain new system-theoretic derivations of the 
Beurling-Lax representation theorem for shift invariant subspaces in both the non- 
commutative and commutative settings; the result for the noncommutative setting 
is due originally to Popescu 03] and for the commutative setting to Arveson jj] 
and McCullough- Trent [^Hj • We also indicate connections with dilation theory and 
the von Neumann inequality for these settings (see |4~5"1 14"81 12~81 17] ) . 

Closely related to the kernels Kc,a(z,w) and K^, A (X,C) (given by l|l.l(j|l and 
(I1.20JI with H normalized to be the identity operator) are kernels of de Branges- 
Rovnyak type (see |3JJ for the classical case) 



(where z — (z%, . . . , z d ) and w — (w%, . . . , Wd) are two sets of noncommuting indeter- 
minates with fcg z (z, w) — ^2 ve:Fd z v w v equal to the noncommutative Szego kernel 
while A = (Ai, . . . , A<j) and C = (Ci, • • • , Cd) are two sets of commuting variables) 
for respective reproducing kernel Hilbert spaces H(Ks), W(iTf) in the respective 
noncommutative and commutative settings. In this situation (where Kg and Kg 
are positive kernels in noncommuting and commuting variable, respectively), the 



K s {z,w) = k Sz (z 7 w) - S(z)k Sz (z,w)S(w)* , 



i-s(\)s(cr 

1-(A,C> 
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respective power series 

S(z) = S v z v , S(X) = J2 S ^ n , 

are contractive multipliers, i.e., the respective multiplication operators 

M s : /(*) -> S(z) ■ f(z), M s : /(A) - S(X) ■ /(A) 

are bounded from H^(!F d ) into Hy(F d ) and from Hu{k d ) into TCy(kd) respectively 
with norm at most 1. A particular issue is the construction of operators 

By.U^X for j = l,...,d and D:U^y 

for some input space IA so that 

S{z) = D + C{I- z x A x z d A d )- x (z x B x + ■■■ + z d B d ), 

S{\) = D + C{I- X X A X \ c iAd)~ 1 (\ x B x + ■■■ + X d B d ) 

satisfy 

%(^) = K s (z,w), K% iA (\,C) = X|(A,C). 
With the resolution of this issue, then the results here lead directly to representa- 
tions of backward-shift-invariant subspaces as reproducing kernel Hilbert spaces of 
the form TL{Ks) and Ti.(Kg) for a Schur multiplier S in both the noncommutative 
and commutative settings as well as linear-fractional realizations for Beurling-Lax 
representers of shift-invariant subspaces for both the noncommutative (see |44p 
and commutative (see |38| ) settings. We work out these issues for the commutative 
setting and for the noncommutative setting in j2] and |10| respectively. 

The paper is organized as follows. After the present Introduction, Section [21 
focuses on the noncommutative Fock space setting while Section [3] focuses on the 
Arveson-space setting. Section [21 is divided into Section |2~T1 dealing with the con- 
nections between solutions of generalized Stein equations and strong stability of 
the state dynamics for noncommutative Fornasini-Marchesini systems and Section 
12. 21 dealing with characterizing ranges of observability operators as backward-shift- 
invariant subspaces of the Fock space with a certain reproducing-kernel-Hilbert- 
space structure. The first subsection f Section 13. 1|1 of Section [3 deals with the less 
tractable issues parallel to the material in Section [2.1l of generalized Stein equations 
and stability for commutative Fornasini-Marchesini systems and also presents the 
abelianization map giving the connection between noncommutative and commuta- 
tive Fornasini-Marchesini systems. The second subsection fSection l3.2[) of Section^ 
parallel to Section 12.21 discusses characterizations of observability-operator ranges 
for the case of a commutative Fornasini-Marchesini state-output system. The re- 
sults are the most satisfying in case the operator-tuple A giving the state dynamics 
is commutative — these are collected in Subsection 13 . 2 1 ll The more implicit results 
for the case of noncommutative A are given in Subsection 13 . 2 . 2l 

2. The Fock-space setting 

2.1. Output stability and Stein equations: the noncommutative case. For 

d a positive integer, let T d be the free semigroup T d generated by the set of d letters 
{1, . . . , d}. Elements of T d are words of the form ijv ■ ■ -i x where ig S {1, . . . , d} 
for each £ = 1 , . . . , N with multiplication given by concatenation. We also use 
to denote the empty word; this serves as the unit element for T d . For v = 
iNiN-i ■■■ii G J~d, we let |i>| denote the number N of letters in v and we let 
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v T := i\ ■ ■ ■ ijv-ii/v denote the transpose of v. We let z = (z\, . . . , Zd) to be a 
collection of d formal noncommuting variables and let y((z)} denote the set of 
formal noncommutative series X^e-F f vZV wnere fv^y and where 

z v = z lN z lN _ 1 ■ ■ ■ z h if v = i N i N -i (2.1) 
The Fock space iy{Fd) is defined as 

:= I {/„ W„ : J|/4<oo L (2.2) 

I v&F d ) 

If we let Xv be the characteristic function of the word v, so 

1 if v' = v, 



Xv = {Xv (v')}v>er d where 



otherwise, 



and we let By be an orthonormal basis for y, then {xvVi ■ v G J-'diUi £ £>;y} is an 
^(J" d ). The space 4( 



orthonormal basis for iy(J-d)- The space ^(-^"rf) can be identified as the tensor 



product l 2 (J-d) <8> y and is mapped unitarily onto the space 

^(^) = {e f^ v ey{{z)): 52\\M\y«x\ (2.3) 
by the noncommutative Z-transform 

{/.w.-n^E// (2.4) 

with the monomials z 11 playing the role of the basis vectors \v ■ 

The noncommutative multidimensional analogue of the system (|1 . 1|> is the system 
with evolution along the free semigroup Td given by {EE} . Upon running the 
system (|1.13|) with the zero input string u(v) = for v £ Td and a fixed initial 
condition x(0) = x 6 X we get 

y(v) = CA v x + CA V " B 3 u(v'). (2.5) 

Here we extend the noncommutative functional calculus (|2.1|l from noncommuting 
indeterminates z = (zi, . . . ,Zd) to a d-tuple of operators A = (Ai, . . . , Ad); we use 
the notation 

A v = A iN A lN _ l ■ ■ -A h if v = i N iN-i ■ ■ -h e T d 

where the multiplication is now operator composition. Application of the formal 
noncommutative Z-transform 12.4fl then gives 

y{z) = C(I - Z(z)A)- 1 x{9) + Tv{z)u{z) (2.6) 

where the formal power series Ts (z) (by definition equal to the transfer function of 
the system H1.13|l ') is given by 

T s (z) = D + C{I- Z{z)Ay 1 Z{z)B 

where we have set 



(2.7) 





'A{ 




~B{ 


Z(z) =[z x ■■■ z d ] ®I X , A = 




, B = 
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For details see [16] or, for a more general setting of structured noncommutative 
multidimensional systems, see 

In analogy to the classical case, the system (I1.13f) is called output-stable (and 
in this case we will say that the pair (C, A) is output-stable) if the output string 
{y(v)}veF d , defined as in (|2.5|l but with the input string {u(v)} v ^jr d assumed to be 
equal to 0, belongs to £y(J-d) for every x 6 X and the observability operator 

O c ,a- x~{CA v x} vf -r d (2-8) 

is bounded as an operator from X into iy{Td)- The Z-transformed version of 0c, A 
is 

Oc,a : x h- (CA v x) z v e y((z)) 

vEFd 

and the following realization formula for Oc.a is immediate: 

Oc,ax = C(I - Z(z)A)- 1 x. 

If (C, A) is output-stable, then Oc,a maps X into Hy(!Fd) and is bounded. In this 
case it makes sense to introduce the observability gramian 

Qc,a := (OcaTO^a = (Oc,a)*Oc,a (2.9) 
and its representation in terms of strongly converging series 

Gc,a = ]T A* vT C*CA v (2.10) 

follows immediately by definition (|2.8|l of Oc,a and the formula (|2.2|) for the norm 
in £y(J-d)- The second equality in l|2.9|) follows by definition of Oc.a and the 
formula (|2.3[) for the norm in Hy{Td)- 

Definition 2.1. A pair (C, A) is called observable if Gc,A is positive definite and 
exactly observable if Qc,A is strictly positive definite. We say that the (i-tuple of 
operators A — (A±, . . . , Ad) is strongly stable if 

lim V ||A"a;|| 2 -> for all x 6 X. (2.11) 

1)6^: |u|=Af 

We mention that the term pure rather than strongly stable has been used in this 
context (see |S]), but we prefer the present terminology since pure so as to avoid 
confusion with the use of the term pure in the context of contractive operator- valued 
functions (see |4L)| ) . 

In analogy with the classical case one can introduce the unobservable subspace 

Ker£c,A = Ker£> c ,A= f] KerCA v . (2.12) 

ve^d 

Thus, observability of (C, A) means that Ker Qc,a is the zero subspace or that 

CA v x = (Vt> e T d ) => x = 0. (2.13) 
The following is the noncommutative Fock-space counterpart to Theorem ll.il 
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Theorem 2.2. Let A = {Ai,...,A d ) £ C(X) d and let C e C(X,y). Then the 
pair (C, A) is output-stable if and only if the (generalized) Stein inequality 

H - A*HA\ A* d HA d >C*C (2.14) 

has a positive semidefinite solution H G C(X). In this case, 

(1) The observability gramian Gc,A satisfies the generalized Stein equation 

H — A\HAx A* d HA d = C*C (2.15) 

and is the minimal positive semidefinite solution of the generalized Stein 
inequality (|2.14|l . 

(2) The positive semidefinite solution of the Stein equation (|2.15|l is unique if 
A is strongly stable, i.e., ()2.11JI holds. Moreover, in case A is contractive 
in the sense that 

A\Ai + >-- + A* d Ai<I Xl (2.16) 

then the solution of the Stein equation (|2.15|l is unique if and only if A is 
strongly stable. 

Proof. Suppose first that (C, A) is output-stable. Then for each x € X, 
{CA v } ve r d e £ 2 y (T d ), i.e., Y \\CA v xf y < oo. 

This has the consequence that the infinite series 

oo 

Yj a v *c*ca v 

converges in the strong operator topology to an operator H 6 C(X) (in fact, H = 
Gc,a is the observability gramian). From this infinite-series representation for Gc,A 
it is easily verified that Gc.a is positive semidefinite and satisfies the Stein equation 
(I2.15|l and hence also the Stein inequality l|2.14|l . 

Conversely, suppose that the Stein inequality l|2.15|l has a positive semidefinite 
solution H. We first claim that 

H> J2 A*" T C*CA V + A* vT HA v (2.17) 

ver d : \v\<N v£F d : |«|=JV+1 

for each feZ + . For N — 0, l|2.17|l collapses to l|2.14|l which is given. Inductively 
assume that 

H> Y A* vT C*CA v + A* vT HA v . 

v&r d : \v\<N u€Jf d : \v\=N 

Use the Stein inequality (|2.14|) to replace H on the right side by its lower bound 
C*C + A\HA\ H h A* d HA d to get from this 

H> Y A* vT C*CA v + Y A* vT HA v + Y A* yT C*CA V 

veT d : \v\<N v£F d : |t,|=V+l v£F d : \v\ = N 

which then simplifies to (|2.17|) as wanted. 
We rewrite (I2.17|) in the form 

Y A* yT C*CA V <H— Y A* yT HA V < H. (2.18) 

v£F d : \v\<N veF d : \v\=N 
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By letting N — ► oo in l|2.18(l we conclude that the left hand side sum converges 
(weakly and therefore, since all the terms are positive semidcfinitc, strongly) to a 
bounded positive semidefinite operator. By 12.11) [I . 

lim V A* yT C*CA V = V A* yT C*CA V = Gc a 



N — * 



DC 



v£F d : \v\<N v£F d 



and passing to the limit in (|2.18|) as N — > oo gives Gc,a < H. In particular the 
operator Gc,a is bounded (since H is) and therefore the pair (C, A) is output- 
stable. □ 

Proof of (1): As observed in the proof of the first part of the theorem, from the 
infinite-series representation (|2.10(l it follows that Gc,A satisfies the Stein equation 
1)2.15(1 . If H is any solution of the Stein inequality, the computation leading to 
(I2.18|l shows that H satisfies (|2.18|) . By taking the limit as N — > oo we conclude 
that Gc,a < H as asserted. □ 

Proof of (2): Suppose that A is strongly stable and that H solves the Stein equation 
(|2.15|l . Then the proof of (|2.17(l shows that in this case 1)2.17(1 holds with equality: 

H= A*" T C*CA V + ^ A* yT HA" (2.19) 

v£F d : \v\<N v£F d : \v\=N+l 

for each N = 0,1,2,.... Taking the limit as N — > oo and using the stability 
assumption 12.11(1 we conclude that H = Gc,a- 

For the converse direction we assume in addition that A is contractive (i.e., 
d2.1tj [I holds). We prove the contrapositive: if A does not satisfy the stability con- 
dition 12.110 . then the solution of the Stein equation 12.15fl is not unique. Assume 
therefore that A is not stable. By the assumption 12.16|) . the sequence of operators 



A N = Yl A**'A*, N = l,2,... 

vef d ■ \v\=N 

is decreasing and bounded below and therefore has a strong limit A. Since A is 
assumed not to be stable, this limit A is not zero. However it is easily verified that 

A\A N A 1 + --- + A* d A N A d = A N+1 . (2.20) 

Taking limits in 12.20fl gives that A = limN^oo An satisfies the homogeneous Stein 
equation 

A - A^AAi A* d AA d = 0. 

We conclude that the solution of the Stein equation (|2.15|l cannot be unique. □ 

Particular cases of output pairs (C, A) are the cases where (C, A) is contractive 
(i.e., the Stein inequality (|2.14|) holds with H = Ix) and where (C, A) is isometric 
(i.e., the Stein equality (|2.15|) holds with H — Ix)- For these cases some additional 
observations can be made along the lines of Theorem 12. 21 

Proposition 2.3. (1) Suppose that (C, A) is a contractive pair. Then (C, A) 
is output-stable with Gc.A < Ix cind the observability gramian Gc.A is the 
unique positive semidefinite solution of the Stein equation 12.15fl if and only 
if A is strongly stable. 
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(2) Suppose that (C, A) is an isometric pair. Then (C, A) is output- stable. 
Moreover H = I is the unique solution of the Stein equation Ij2.15|l if and 
only if A is strongly stable. In this case Oc,a is isometric and hence also 
(C, A) is exactly observable. 

Proof. Statement (1) immediately follows from statements in Theorem 12.21 com- 
bined with the observation that (C, A) being a contractive pair implies that A is 
contractive. 

The first two assertions in statement (2) follow in a similar way. As for the 
last assertion, for the case where (C, A) is isometric, lx is a solution of the Stein 
equation (|2.15(l : for the situation where A is strongly stable, uniqueness implies 
that the observability gramian Gc,a — Ix, i-e., that 0c, A is isometric. Then also 
(C, A) is exactly observable by definition. □ 

Remark 2.4. The converse of the last part of Proposition 12 . 31 does not hold even 
for the case d = 1 . More precisely, there exists an isometric pair of operators (C, A) 
such that (C, A) is observable but A is not strongly stable. 

An example necessarily requires that dim X = oo. In the terminology of Sz.- 
Nagy-Foia§ |4(J) . it suffices to produce a completely non-isometric (c.n.i.) contrac- 
tion operator A on a nontrivial Hilbert space X (so dim X > and there is no 
nonzero- invariant subspace M for A such that A\m is an isometry) in the class 
C\. (so A n x — » in X for some x € X implies that x = 0). Indeed, if A is 
such an operator, set C = (I — A* A) 1 / 2 considered as an operator from X into 
y := Ran(J - A* A) 1 / 2 (the closure of the range of (I - A* A) 1 / 2 ). Such an A is not 
strongly stable by the definition of the class C\., the definition of C makes the pair 
(C, A) isometric, and the condition that A is c.n.i. implies that (C, A) is observable. 

To construct such an operator A, let 9 be a Schur-class outer function such that 
log(l — \6\ 2 ) is not integrable (with respect to arc- length Lebesgue measure) over 
the unit circle T. Furthermore, let K(0) be the associated Sz.-Nagy-Foia§ model 
space 



K(0) 



H 2 {B) 
L 2 (T) 



e 



0(A) 

(i- 10(C)! 2 ) 1 / 2 



H 2 (B) where A e D and ( e T, 



and let S(0) be the Sz.-Nagy-Foia§ model operator 



T — Pk(6) 



M x 
M c 



K(0) 

where M\ and are the operators of multiplication by A and by C, respectively. 
Now we let A := S(0)* and note that A is in the class C\. by Proposition 3.5 
in 00] (since 9 is outer) and A is c.n.i. by Theorem 5 in J^j (since the non-log- 
integrability property of I — \0\ 2 implies that there is no ff°°-function a(z) for which 
|a(C)| 2 < 1 - |0(C) I 2 for C S T). This completes the construction. 

Let us say that the pair (C, A) is similar to the pair (C, A) if there is an 
invertible operator S on X so that 

C = CS~\ Aj = SAjS- 1 for j = 1, . . . , d. 

Then we have the following characterization of pairs (C, A) which are similar to a 
contractive or to an isometric pair. 
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Proposition 2.5. (1) The pair (C, A) is similar to a contractive pair (C, A) 
if and only if there exists a bounded, strictly positive-definite solution H to 
the Stein inequality <|2.14[1 . 
(2) The pair (C, A) is similar to an isometric pair if and only if there exists a 
bounded, strictly positive-definite solution H of the Stein equation (|2.15|l . 

Proof. Suppose that H is a strictly positive-definite solution of H2.14[l . Factor H 
as H — S* S with S invertible and set 

C = CS-\ Ak^SAkS- 1 for k = l,...,d. (2.21) 

Multiplying l|2.14|l on the left by and on the right by 5 _1 then leads us to 

I -AIM A* d A d >C*C, 

i.e., (C, A) is a contractive pair which is similar to the original pair (C, A). Con- 
versely, if (C, A) given by l|2.21|l is contractive, then H = S*S is bounded and 
positive-definite and satisfies the Stein inequality H2.14JI . This verifies the first 
statement of the Proposition. The second statement follows in a similar way. □ 

As a consequence of the observations in Proposition ^. 51 Proposition 12 . '6\ can be 
formulated more generally as follows. 

Proposition 2.6. (1) If the pair (C, A) is such that the Stein inequality (|2.14(1 
has a strictly positive-definite solution H , then (C, A) is output- stable. 
Moreover, the observability gramian Gc,A is the unique positive semidefinite 
solution of the Stein equation (|2.15|1 if and only if A is strongly stable. 
(2) // the pair (C, A) is such that the Stein equation (I2.15|) has a strictly 
positive- definite solution H , then (C, A) is output-stable and the observabil- 
ity gramian Qc,A is the unique positive semidefinite solution of the Stein 
equation H2.15JI if and only if A is strongly stable. In this case (C, A) is 
moreover exactly observable. 

The last part of Proposition 12.61 has a converse. 

Proposition 2.7. Suppose that the pair (C, A) is output-stable and exactly observ- 
able. Then A is strongly stable, i.e., (|2.11(l holds. 

Proof. If (C, A) is output-stable and exactly observable, then the observability 
gramian Gc,a is a strictly positive-definite solution of the Stein equation (|2.15|l . 
Hence (|2~T§j) holds with H = g c ,A- 

(Qc,ax,x)= {A.* vT C*CA v x,x)+ (Qc,aA v x,A v x). 

veF d : \v\<N uefj: M=JV+1 

(2.22) 

From the infinite-series representation (|2.10|) for Qc.a, taking limits in (|2.22(l gives 
lim V (g c ,AA v x, A v x) = 0. (2.23) 

ver d : \v\=N+l 

The strict positive-definiteness of Qc.a tells us that there is an e > so that 

e\\x\\ 2 < {Qc,ax,x) for all x e X. (2.24) 
In particular, from (|2.24l) with A v x in place of x combined with 12.23f) we get 

e J2 II A ^H 2 < E (Gc,aA*x,A v x)^0 

v£F d : |u|=iV+l v£F d : \v\=N+l 
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for all x G X, and we conclude that A is strongly stable as asserted. □ 

2.2. Observability-operator range spaces and reproducing kernel Hilbert 
spaces: the noncommutative-variable case. To develop the noncommutative 
analogue of Theorem 1 1.21 we first introduce the right noncommutative shift opera- 
tors Sf , . . . , Sf on H^(T d ) as follows: 

Sf : f" zV " E f« zV *i = E fv* Vj C? = !» • • • > d Y ( 2 - 2 5) 
vefd veF d ve^d 

It is readily seen that their adjoints (backward shifts) are given by 

(Sf)* : £ ' * E ^ C? = 1, • • • > d). (2.26) 

Their left counterparts Si,...,S%, also on Hy{!Fd), are given by 

E f- z " " E M'*" = E (2.27) 

uGJ^d ^G^d ■u6^ r d 

with adjoints given by 

{sf)*-- E /^"^ E fs** v - ( 2 - 28 ) 

Let t denote the unitary involution on Hy(J-d) given by 

r: Y, ^ - E /» T *°- (2-29) 
»efd «e^d 

In addition to the unitary property r* = r _1 of r, note also that r intertwines the 
left shifts with the right shifts: 

{Sf)*T = r(Sf)*, Sfr = rSf for j = 1, . . . , d. (2.30) 

Then we have the following Fock-space analogue of Theorem 1 1.21 

Theorem 2.8. Suppose that (C, A) is an output-stable pair. Then: 

(1) The intertwining relation 

(Sf)*d c ,AX = d c ,AAjX (xeX) (2.31) 

holds for every backward- shift operator (Sj)* defined in (|2.26|) and hence 

Ran Oc,a is (S* 1 )* -invariant for j — 1, . . . , d. 

(2) Let H > be a solution of the Stein inequality l|2.14|l and let X' be the 
completion of X with H -inner product \\[x\\\ 2 x , = (Hx,x)x- Then Aj and 
C extend to define bounded operators Aj : X' — > X' for j = 1, . . . , d and 

C : X' — > y and the observability operator Oc,a extends to define a con- 
traction operator Oc'.A' from X' into Hy{Td). Moreover, Oc>,A>'- X' — > 
Hy(J-d) is an isometry if and only if H satisfies the Stein equation (|2.15(l 
and A' = (A[, . . . , A' d ) is strongly stable, i.e., 

Y {HA v x,A v x)^0 for all x eX. (2.32) 

v£F d : \v\=N 
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(3) If H > is a solution of the Stein inequality (|2.14(l and the linear manifold 
M. := RanOc.A is given the lifted norm 

\\d c ,Ax\\ 2 M = „inf _ (Hy,y) x , (2.33) 

yeX: Oc,aV=O c ,ax 

then 

(a) M. can be completed to M.' = RanOc,A' (with (C',A') as in #2 
above) with contractive inclusion in Hy(Fd): 

\\.f\\ 2 H U r d) < \\f\\ 2 M> for all f g M'. 

Furthermore, M! is isometrically equal to the formal noncommutative 
reproducing kernel Hilbert space with reproducing kernel Kc,a-h given 
by HHBJ>. 

(b) The following difference- quotient inequality is valid 

E WiSfT f\\ 2 n { K c ^ H) < \\f\\ 2 H(K c , Aia) ~ Whfy (2-34) 

3=1 

for every f G M! = TI{Kc,A;h) with equality holding in 12.34JI if and 
only if H2.15(l holds. 

(4) Conversely, if AA is a Hilbert space included in Hy(Td) which is invariant 
under (S^)* for j = 1, . . . , d such that the difference- quotient inequality 

d 

£ll(sf)7lft<< II/ILm-II/0]]* (2.35) 

3=1 

holds for every f g M., then M is contractively included in Hy(Td) and 
there exists a contractive pair (C, A) (so H = I positive definite solution 
of the Stein inequality i|2.14[l ) such that 

M = H(K c ,A;i) =RanO c , A 

isometrically. In case (|2.35() holds with equality, then (C, A) can be taken to 
be an isometric pair. An explicit (C, A) meeting these conditions is given 
as follows. Take X to be the Hilbert space X = t(M) (where r is the 
involution given by (|2.29|) ) with ||t(/)||;k = ||/||ai and define C: X — > y 
and A = (Ax, . . . , Ad) on X by 

A 3 = (Sf)*\ x forj = l,...,d; C = E\ x :X^y, H = I x (2.36) 
where E: Hy(J-d) — > y is given by 

E: fvz v ^h- (2.37) 

Proof of (1): Applying [Sf 1 )* to a typical element from Ka,nOc,A (with notation 
as in H2.7JI ). we get 

(S*r(C(I-Z(z)A)-i X ) = (Sf)*(£(CA^) 

= {CA. vl x)z v 

= C(I - Z(z) A)- 1 A jX eR&ndcA- (2.38) 
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The latter equality shows that Ran Oc,a is invariant under (Sj)* for all j — 1, . . . , d 
(backward-shift-invariant) and Ij2.31|) follows. □ 

Proof of (2): The Stein inequality H2.14[) amounts to the statement that (C, A) 
is contractive and well-defined on the dense subset [X] of X' (where [x] is the 
equivalence class containing x) and hence extends to a contractive pair (C",A') 
on all of X' and moreover the inequality (|2.17(l holds for all N — 1,2,.... From 
this we see that Oc,a is contractive from X with the 7J-inner product to Hy{Td), 

and hence also 0c, A' is contractive from X' to Hy(Td)- The inequality H2.17JI is 
actually a chain of inequalities Wm > Wm+i for N = 1,2,..., where 

W N = A* vT C*CA v + A* yT HA V . 

vE^d- M<JV-1 u€Jf d : |«|=JV 

Note that 

s-limjv^oo Wn = Qc,a + Ah,a 

where 

A^.A = s-limAf^oo ^2 A * V HA - 

\v\=N 

In particular it follows from l|2.16|) that 

d 

H > C*C + AjHAj > W N for all N > 2 
i=i 

and hence, by taking the strong limit on the right hand side, we get 

d 

H> C*C + Y,A*HA,j > £c,a + A h , a . (2.39) 

i=i 

By definition, Oc'.A' ■ X' — > Hy^d) being an isometry means that Gc.a = H in 
which case (|2.39ll becomes 

d 

H > C*C + A*jHAj > H + Ah,a (2.40) 
i=i 

which in turn forces Ah. a = and equalities throughout (|2.4(J|) . The condition 
A^a = just means that A' is strongly stable. From equality holding in (|2.4(J|) 
we see that the Stein inequality (|2.14|l holds with equality, i.e., the Stein equation 
(I2.15|l holds. Conversely, by reversing the steps of the argument, we see that A' 
being strongly stable and the Stein equality holding leads to Gc,a — H, i.e., to 
Oca' being an isometry from X' into Hy{Td)- □ 

Proof of (3a): Statement (3a) follows from general principles laid out in ^5]. □ 

Proof of (3b): For / of the form f(z) — C(I — Z(z)A)~ 1 x, we have 

Wf\\n(K C ,A:H) 



< 2 - (Hx,x). 



We see that then 

|2 



S VWh{k 0:A . iH ) = (HA 3 x, Ajx) for j = 1, . . . , d (from (239), h - Cx. 



18 J. A. BALL, V. BOLOTNIKOV, AND Q. FANG 

With these substitutions, we see that 1(2. 34J) is equivalent to 



J2(HA jX ,A 3 x) x < {Hx,x) x - \\Cx\\ 

3=1 

or, in operator-theoretic form, 

AlHAi + ■■■ + A* d HAd <H- C*C 



2 



(2.41) 



with equality in 1(2.34(1 equivalent to equality in 1(2 .41(1 , This completes the verifi- 
cation of part (3b) of Theorem 12. 81 □ 

Before commencing the proof of part (4) of Theorem 12. 81 we collect some useful 
facts concerning Hy(J-d) itself. 

Proposition 2.9. Let S = (Si, . . . , Sd) denote either the right shift S R or the left 
shift S L 



(Si , ■ .., Sd), 



defined as in (|2.25|l and (|2.27ll and let the operator E: Hy(Td), — > y be defined 



in (|2~37|) . Then: 

(1) The operator-tuple S* = (S*, . . . , S%) is strongly stable, i.e., 

lim V ||S^/l| 2 = for each f G H y (T d ). 
ver d -. \v\=n 

(2) The operator 



N- 



(2.42) 



Hl(F d ) -» {Hy(T d )) d ®y 



is unitary, i.e., 

EE*=I y , ES % =0, S[S, 6, J for i,j = l, 
(where 5ij stands for the Kronecker symbol), and 
I — S\Si — • • • — SdSd — E E. 
(3) X — I is the unique solution of the Stein equation 
X — S\XSi — ■■• — SdXSd — E* E. 



(2.43) 
(2.44) 
(2.45) 



(4) For every f G H y (T d ), 



f(z)-h=J2(SjS*f)(z 

3=1 



3=1 
d 



(2.46) 



• ((Sffh)(z) i fS = S> 

(3=1 



(5) The observability operator O e ga. is equal to the operator r defined in 
(12.29(1 and hence is unitary. 
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Proof. If /(2)=£„ e ^/^,then 

(S R * v 'f)(z) = ^ (S w /)(*) - E /f' T **" 

and hence, in either the left or the right case, we have 

ES* v f = f v T. (2.47) 

Therefore, 

E H S *VI| 2 = E UMI 2 ^ aS ^-°°> 
uS^: |«|=JV uS^: \v\>N 

and follows. Equalities (j2~4l|) and (|2~4l|) follow from (|2~2"5jl . lf2~2"H|l . l(2~27jl . 

(|2.28fl l|2.37f) and the fact that £7* is the inclusion map of y into H y {Td)- Applying 
the operator identity l|2.44() to an / e H y (T d ), we get (|2.46() . Finally, from (|2.47|) 
we see that, for both the left and the right case, 

d E , S 'f = E (#S*VK - E ^ TZ ' = T ^ 

for all / G Hy^d)- That X = / is the unique solution of the Stein equation (|2.45|) 
is now a consequence of (|2.42|l combined with the last part of Theorem 12. 21 □ 

Proof of (4) in Theorem \2.8t Suppose that M. is a Hilbert space contractively in- 
cluded in HyiTd) which is invariant under (S^)* for each j = 1, . . . , d such that 
the difference-quotient inequality l|2.34[l holds. Set X = r(Ai) with norm inher- 
ited from A4. From the intertwining relations (|2.30() we see that X is invariant 
under the left backward shifts (5f)*, . . . , (5^)*). Define operators A j : X — > X 
for j = 1, . . . ,d and C : X — > 3^ by (|2.36(l . From the difference-quotient inequality 
(|2.35|) together with the definition of the ,Y-norm and the intertwining relations 
(|2.3U[) , we have 

d d d 

E \\( s tr(Tf)\\% = E \wsjTf)\\x = E \\(s?rf\\ 2 M 

i=i j'=i i=i 

<ll/lll,-ll/0ll^ = H/)lli-ll(T/)0ll^ 

and hence H — Ix satisfies the Stein inequality (|2.14|) . From Proposition 12.91 we 
see that 

Cc,a = O e ^l,\x = t\x= t (M) 

and hence Oc.at\m = Im- Therefore, for each / £ M. we have 

ll/lk(K C , A;l) = 0c,at/||„ ( k c , a ) = l|r/|U = II/IIm 

and thus Ai = H(Kc,a) isometrically. It then follows from part (3a) of the theorem 
that in fact M. is contractively included in Hy(Td)- □ 

As explained by part (4) of Theorem 12. 81 for purposes of study of contractively- 
included, backward-shift-invariant subspaces of Hy^d) which satisfy the difference- 
quotient-inequality H2.34f> . without loss of generality we may suppose at the start 
that we are working with X' as the original state space X and with the solution H 
of the Stein inequality (|2.14(1 to be normalized to H = Ix- Then certain simplifi- 
cations occur in parts (l)-(4) of Theorem 12. 81 as explained in the next result. 
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Theorem 2.10. Suppose that (C, A) is a contractive pair with state space X and 
output space y . Then: 

(1) (C, A) is output-stable and the intertwining relation (|2.31|l holds. Hence 
Ran Oca is invariant under the backward shifts (S^)* for j = 1, . . . , d. 

(2) The observability operator Oq.a is a contraction from X into Hy(J-d). 
Moreover Oc,a is isometric if and only if (C, A) is an isometric pair and 
A is strongly stable. 

(3) If the linear manifold M := Ran Op a is given the lifted norm 

\\dc,Ax\\ M = \\Qx\\x (2.48) 

where Q is the orthogonal projection of X onto (Ker Oc,a) > then Oc,a 
is a coisometry of X onto M. Moreover, A4 is contained contractively in 
Hy(J-d) and is isometrically equal to the formal noncommutative reproduc- 
ing kernel Hilbert space TL{Kc,a) with reproducing kernel Kc,a(z, w) given 
by 

K c ,a(z,w) = C ( J - Z(z)A)~ 1 (I - A*Z(w)*)- 1 C*. 

(4) If Oc,a is given the lifted norm || • ||h(-R" c a) as * n Q2.48|L then the difference- 
quotient inequality 

E \\(S*rf\\ 2 H (KG , A) < \\n 2 H {Kc , A) WUfy (2-49) 

holds for all f G H(Kc.a)- Moreover, (I2.49|) holds with equality if and 
only the orthogonal projection Q of X onto (Ker O^a) 1 satisfies the Stein 
equation 

d 

Q-J2 A *jQ A i = C * C - (2-50) 

In particular, if (C, A) is observable, then (|2.49|) holds with equality if and 
only if (C, A) is an isometric pair. 

Proof. Statements (l)-(3) and all but the last part of statement (4) are direct 
specializations to the case H = Ix of the corresponding results in Theorem 12.81 It 
remains only to analyze the conditions for equality in (|2.49(l . 

From the intertwining relation Ij2.31|l . we see that equality in H2.49J) for a generic 
element / = Oc,a% £ H(Kc,a) means that 

d 

E W d cMjx\\n(K c , A) = \\dc,Ax\\ 2 niKO:A) - ||Cz||5, 

3=1 

for all x G X. Using the definition (|2.48|) of the H(Kc, A)-norm, we rewrite this 
last equality as 

d 

\\n A .^.||2 _ ||n^l|2 _ ll/-7~l|2 



\\QA j x\\ x = \\Qx\\ x -\\Cx\\t 
This holding for all x G X is finally equivalent to the Stein equation (|2.50ll . □ 



3=1 
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Remark 2.11. In Theorems 12.81 and 12 . 1 01 we could equally well have interchanged 
the roles of left versus right. For a given output pair (C, A), define the associated 
left observability operator A : X — > Hy(J-d) by 

O h c>A .x= £ CA»V. 
Then the linear manifold Ran Oq a is invariant under the left backward shifts 



{St 



(Sfi)*) as verified by the intertwining relation 



We leave the precise statements and proofs to the interested reader. 

The characterization (|2.50() of the difference-quotient inequality holding with 
equality for a space TL{Kc,a) in Thcor cm !2 . 1 01 can be made more explicit as follows. 

Proposition 2.12. Suppose that (C, A) is a contractive pair as in Theorem \2.1(A 

and let Q be the orthogonal projection onto (Ker Oc,a) ■ Then Q satisfies the Stein 
inequality 

Q-A* 1 QA 1 A* d QA d >C*C (2.51) 

and we have the inequalities 

Gc,A<Q<Ix. (2.52) 
If we write C,Aj,Q in 2 x 2-block matrix form with respect to the decomposition 
X = KerO CiA © (Ker (5c, a)- 1 as 



C = [0 C*°] 



'A n A j2 ~ 




"o 


o" 


. A°_ 


Q = 





i 



(2.53) 



for j = X,...,d, then Q satisfies the Stein equation (|2.50() if and only if the 
pair (C , A ) is an isometric pair, in which case we also have that Aj2 — (so 

(Ker Oca) 1 is invariant for Aj) for j = 1, . . . , d. 

Proof. First note that Ker Oc,a is invariant for each Aj and that Ker Oc, a C 
Ker C. Therefore the matrix decompositions of C, Aj , Q with respect to the de- 
composition X — KerOc^A ® (Ker Oca)" 1 have the form as given in (|2.53(l . Next 
note that the contractive property of the pair (C, A) means that 





C°*C° 



d 

£ 

3=1 



A* x A 3l A* jX A j2 
A* 2 A jX A* j2 A j2 + AfAf 



< 


I 


0" 












On the other hand, the Stein inequality (|2.51(l works out to be 

d r 





C°*C° 



E 

3=1 



A* j2 A n AfAf 



< 



I 
I 



(2.54) 



(2.55) 



As the left hand side of (|2.55|) is dominated by the left hand side of l|2.54|l . it is 
clear that (|2.55|) follows from l|2.54|l , and hence (|2.51|) holds as asserted. Since Gc,A 
is the minimal positive semidefinite solution of the Stein inequality l|2.14|l (by part 
(2) of Theorem 12 .2(1 and we now know that Q is one such solution, it follows that 
Gc,a < Q- As Q is an orthogonal projection on X, we also have Q < Ix and (|2.52|l 
now follows. 
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From the (2,2) entry of Ij2.54jl . we read off 



+ E A h A v + E A ? A ° < Wso, A) ^ 

3=1 3=1 



(2.56) 



In particular 



(KorO c , A )-L' 



i.e., Q satisfies (12.511) . On the other hand, the validity of (|2.5U|) reduces to 

d r,, 

E 

3=1 





c°*c° 









I 



or simply to 

c *c + E<4 = V r a c , A) - (2-57) 

3 = 1 

Thus, the validity of 12.50|) is equivalent to (C°,A°) being an isometric pair, in 
which case we also have that Aj2 = 0. □ 

Finally, we have the following uniqueness result. 

Theorem 2.13. Suppose that (C, A) and (C, A) are two output-stable, observable 
pairs realizing the same positive kernel 

K c ,a(z,w) := C(I - Z(z)A)-\l - A*Z{w)*)~ 1 C* (2.58) 

= C(I - Ziz)!)- 1 ^ - A*Z{w)*)- 1 C* =: K 5 A (z, w). 

Then (C, A) and (C,A) are unitarily equivalent, i.e., there is a unitary operator 
U : X — > X such that 

C = CU and Aj^U^AjU for j = 1, . . . , d. 

Proof. For any two words a, (3 € equating coefficients of z a w 13 in 12.58fl gives 

CA a A* l3 C* = C'A a A* p C*. 
Hence the operator U defined by 

U: A* f3 C*y i-> A* p C*y (2.59) 
extends by linearity and continuity to define an isometry from 
V v = spm{A*PC*y :peF d ,yey} 

onto 

Tlu =span{A* /3 C*y: e T d , y e y}. 

The observability assumption implies that T>u — X and TZjj — X; hence U : X — » X 
is unitary. From (|2.59() it is easily seen that 

UC* = C* and UA* = A*U for j = 1, . . . , d. 

Since U is unitary we then get 

CU = C and AjU = U A 3 for j = 1, . . . , d 

and we conclude that (C, A) and (C, A) are unitarily equivalent as desired. 
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□ 

2.3. Applications of observability operators: the noncommutative set- 
ting. As an application we give a proof of the Beurling-Lax theorem for the Fock- 
space setting originally given by Popescu |44|. We shall in fact prove a more gen- 
eral version of the Beurling-Lax-Halmos theorem for contractively-included (rather 
than isometrically included) subspaces of Hy^d) due in the classical setting to de 
Branges (see [23)- Our proof is similar to that in 01] but highlights more explic- 
itly the role of an associated observability operator. For this purpose we say that 
a formal power series 9(z) = Y^veF ® vZV ^ ^(^> i s a contractive multiplier, 

also written as 9 is in the d- variable, noncommutative Schur-class S nc .d(U,y), if 
the operator Mg of multiplication by 9 

M e : f(z)^9(z)-f(z) 

defines a bounded linear operator from H^(Td) to Hy(J-d) with operator norm at 
most 1. Such a formal power series 9(z) is said to be inner if moreover the operator 
Mg from Hu{Td) to Hy(J-d) is an isometry 1 . 

Theorem 2.14. (1) A Hilbert space M. is such that 

(a) M. is contractively included in Hy(J-d), 

(b) M. is invariant under the right shift operators S^, . . . , S^: 

SfM C M, 

(c) the d-tuple 

S M = i S M,n ••• j S M,d) where S M,j : = s f\M for j = 1, . . . ,d 
is a row contraction 

^M,l(^M,lT + ' ' ' + ^M,d(^M,d)* — I-M, 

and 

(d) (Sjjf)* is strongly stable, i.e. 

lim V ||(S^r/llM^0 for all f E M, 

n — >oo z — ' 

v^Fd ■ \v\=n 

if and only if there is a coefficient Hilbert space IA and a contractive multi- 
plier £ S nc ^d(U,y) so that 

M = 9-H*(F d ) 

with lifted norm 

P-f\\M = \\Qf\\HZm (2-60) 

where Q is the orthogonal projection onto (Ker Mg) 1 - . 
(2) The subspace Ad in part (1) above is isometrically included in Hy^d) if 
and only if the associated contractive multiplier 9 is inner. 

Proof. We first verify sufficiency in statement (1). Suppose that A4 has the form 
KA = 9 • Hy{Td) for a contractive multiplier 9 with A'f-norm given by p.fiOjl . 
From the fact that \\Mg\\ < 1 it is easily verified that \\9 ■ f\\H^(^ d ) < ||# • f\\M, 
i.e., (a) holds. From the intertwining property Sf'Mg — MgSf (note that 



We prefer to define inner to be isometric rather than partially isometric as in 1441 . 
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is multiplication by Zj on the right while Me is multiplication by 9 on the left), 
property (b) follows. 

If Q is the orthogonal projection onto (Ker Mg) 1 - C H^Td), then the intertwin- 
ing property SfMe — MeSf implies that 

QS? = QS?Q and {Sf)*Q = Q{Sf)*Q for j = !,...,<£ 



(2.61) 



Thus 



^M.d] 



fl 
fd 



M 



< 



e[s? 



Q[s? 



Qh 



Qf d 



sS] 



SS] 



Qh 
Qfd. 



M 



2 



e-fd 



M" 



and property (c) follows. Finally, a short computation shows that 

(s& tj y -.e-f^e- s?*Qf, (sft )*": e-f^e- s R * v Qf 



and hence 



E i 

veF d : \v\=n 



/ii 



2 

M 



E w sR * v Qf\\k 

vEFd : \v\=n 







as n — > oo, and property (d) follows as well. Moreover, if 9 is inner and Ai = 
9 ■ Hy(Fd) with the lifted norm (|2.6Q(I . it is clear that Ai is contained in Hy{J-d) 
isometrically. This completes the proof of sufficiency in Theorem 12. 141 

Suppose now that the Hilbert space Ai satisfies conditions (a), (b), (c), (d) in 
statement (1) of Theorem l2.14l Define a d-tuple of operators A = {A\, . . . , Ad) on 
M by 

Aj = (S^)* for j 1 '/• 



where we use hypothesis (b) to set Sj^ := 
coefficient Hilbert space U so that 



S?\m f° r 3 = 1, • • • , dj and choose the 



dimU = rank(7 - A\A X A* d A d ). 

By hypothesis (c) we may then choose the operator C : Ai — » U so that 

C*C = 1- A\A X A* d Ad- 
Then (C, A) is an isometric pair and, by hypothesis (d), A* is strongly stable. Thus 
by part (2) of Proposition 12.31 it follows that the observability operator 

O c ,a: f^(I-Z(z)A)- 1 f 

is an isometry from Ai into H^(JFd). As observed for the general case in part (1) 
of Theorem 12. 81 we have the intertwining condition 



(sfyd c ,A = dcAs R 



M,j) 
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Taking adjoints then gives 

{dc,x)*Sf = S^ j {d c ,A)*. (2.62) 

Let us set 

where t: M — » Hy(Td) is the inclusion map. From hypothesis (a) that ||t|| < 1, we 
see that ||0|| < 1. From the intertwining relation (|2.62[1 (together with hypothesis 
(b)) it follows that 

esf = sfe 

and it follows (see e.g. j3U) that is a multiplication operator, i.e., there is a 
contractive multiplier 9 <G S nc ^d(U,y) so that = M#. From the fact that 
Oca: A4 — > Hy(Fd) is an isometry, it follows that Ran(Oc,A)* = M and also 
that M. = 9 ■ Hy^d) with .M-norm given by (|2.60l) . This completes the proof of 
necessity in statement (1) of Theorem 12 . 141 for the general case. 

We now consider statement (2). In case A4 is isometrically included in Hy^d), 
for any / € H^(Td) we have 

\\(d c ,A.ysff\\ M = \\sf(d c ^rf\\H^ d ) = \\{6c,AYf\\ Hl{ r d) 

for j = 1, . . . , d, since S R is isometric on Hy(J-d). Since, as was observed above, 
Oc,a is isometric, it follows that 

ll P Ra„a c , A ^/H = ll/H fOTa11 /eRanOc.A 

and hence Ran Oca is invariant under Sj for j — 1, . . . ,d. As Ran Oca is also 
invariant under (Sj 1 )* for each j by (|2.62|l . we conclude that RanOc.A is reducing 
for S R . Since RanC is dense in U by construction, we are now able to conclude 
that Ran Oca IS a ^ °f ^u(-^d) and hence Oc,a '■ M — > Hy(J-d) is actually unitary. 
It then follows finally that = t o (Oc,a)* is isometric and hence 9 is inner as 
asserted. This completes the proof of Theorem 12. 141 □ 

A second application of these ideas is to operator model theory. For this appli- 
cation we are given only an operator-tuple T = (Tj., . . . ,Td) € £(H.) d which is a 
row contraction, so I - TiT x * T d T* d > 0. Set 

Dt* := {I - TiTj* T d T* d ) l!2 and y := RanL> T *- (2.63) 

We apply the ideas of the previous sections concerning the general pair (C, A) to a 
pair of the special form (Df ,T*). For simplicity we assume in addition that T* 
is asymptotically stable, i.e., 

lim V IIT*^!! 2 = for all x G H. 

v£F d : \v\=N 

Then we have the following dilation result. 

Theorem 2.15. Suppose that T = (Ti, . . . , Td) is a row contraction with T* asymp- 
totically stable as above and define the defect operator Dt* and the coefficient space 
y as in H2.63J) . Then there is a subspace M C Hy(J-d) invariant for the backward 
shift operator-tuple S R * on Hy(Td) so that T is unitarily equivalent to PmS r \m- 
In particular, T has a row-shift dilation unitarily equivalent to S R on Hy(J-d). 



26 



J. A. BALL, V. BOLOTNIKOV, AND Q. FANG 



Proof. By the same arguments as in the proof of Theorem l2.14l we see that 
is isometric and satisfies the intertwining 
If we then set 

M = Rand DT „ >T ., 

then Od t , ,t* implements the unitary equivalence between T and Pm^>\m as 
wanted. □ 

Remark 2.16. In the classical case d = 1, the procedure for constructing the 
unitary dilation of a contraction operator via the observability operator as in the 
proof of Theorem 12. 151 corresponds to the construction of Douglas (see [23 ) ( see 
also 001 Section 1. 10.1]) which is an alternative to the more popular Schaffer-matrix 
construction of the unitary dilation (see (401 Section 1.5]). Popescu (see (411 I42p 
used an analogue of the Schaffer-matrix construction to construct the row-unitary 
dilation of a row-contraction operator-tuple. From the existence of this dilation, he 
went on to verify a von Neumann inequality (see |45|): 

\\p(Ti,...,T d )\\ < \\p(Si,...,S d )\\ 

for any polynomial p G C(z) in the noncommuting variable z = (zi, . . . ,z d ). He 
returned to this topic in |48| to give another proof of the von Neumann inequality 
(actually a more general version involving nonanalytic polynomials) based on the 
Poisson transform: for T a strict row-contraction (one can reduce the general case 
of a row-contraction to the case of a strict row-contraction via a limiting procedure) , 
one defines the Poisson transform P(T) : Cffly^a), H) — > C(H) by 

P(T)[X] = (d D ^ tT .)*x6 D ^,T*- (2.64) 

It is argued in 0H| (as well as in 122 in the context of the classical case) that 
this is an elementary (i.e., dilation- free) proof of the von Neumann inequality. In- 
deed, as argued in [22 , this proof of the von Neumann inequality goes back to the 
paper of Heinz However, we would argue that the dilation is very near the 

surface in this proof as well, since the Poisson kernel, i.e., the observability opera- 
tor Od t , ,t* , provides the factorization of the Poisson transform ((2.641) and is also 
the operator embedding the state space H into the dilation space Hy{Td) in the 
Douglas approach to dilation theory. 

3. The commutative-variable Arveson-space setting 

3.1. Output stability and Stein equations: the commutative-variable case. 

To introduce the commutative multidimensional counterpart of the Hardy space 
H 2 (H>), we recall standard multivariable notations: for a multi-integer 

n = (m, ...,rij)6Zf 

and a point A = (Ai, . . . , A<j) S C d , we set |n| = m+n2 + - ■ -+nd, n ' = ni\ri2\ . . . rid). 
and 

A n = A' 1 ll A' 2 l2 ...A" £i . (3.1) 
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The space can be derived from the full Fock space by "letting the variables com- 
mute". For this purpose we introduce the abelianization map a: Td — > Z+ given 

by 

a(i N ■ ■ ■ ii) — (ni, . . . ,n d ) where = it = k} for k = l,...,d. 

A key combinatorial fact is that 

#a -i (n) = M. (3.2) 

We then consider the symmetric Fock space ly(STd) equal to the subspace of 
iy{Td) spanned by the elements XnV (n e 7L\ and y e y) where Xn is given by 

Xn = ^ ' Xv- 

v : a(v)— n 

Note that 

IIXn|l| ( ^= E 1 = ^ 

vefd ■ a(«)=n 

and hence, if £> is an orthonormal basis for J 7 , then an orthonormal basis for ^(STd) 
is the set 

It is then natural to identify t y {STd) with the weighted sequence space l? w yi^\.) 
consisting of all J^-valued Z+-indexed sequences {/n} n ez d f° r which the norm given 

by 

IK/n}nez* ll/i i3 ,(z*) = E ^Wll/nll 2 where w(n) = — , 

neZ d 

is finite. We abbreviate H? w C (Z+) to ^(Z+) and observe that 

Z +)=OZ+)®3>- (3-3) 
The commutative (i-variable Z-transform 

{/n} neZ , ~ r(A) = ]T / n A n 

neZ d 

maps ^(Z+) unitarily onto the Arveson space 




H{k d ) := J /(A) = ]T / n A n : ll/ll 2 = E |St • 1^ 



2 < oo 



with inner product given by 

(f,g)n(k d ) = E l^TT^™^" 



if 

/(A) - E and 5(A) - E 5nA n . 
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Then it follows that the set {^/^A n : n e Z^} is an orthonormal basis for P w (7L d + ). 
By general principles concerning reproducing kernel Hilbert spaces we see that 
7i{kd) is a reproducing kernel Hilbert space of functions analytic on the unit ball 

d 

l d = {A=(Ai,...,A d ) eC d : ^|A fe | 2 <l} 



fc=i 



with reproducing kernel fcd(A, C) given by 



1-<A,C) 



(see e.g. [7]). This justifies the notation Tt(kd) for the space. In analogy to (|3.3[) we 
will use notation Hy(kd) := Tt(kd) <8> y for the tensor product Hilbert space that is 
characterized by 



n y (k d ) = J /(A) = J2 /n A " : Il/H 2 = E ^TT • H/»l 



< oo 



If we define the map n by 

II: {fv}ver d 



E /« 

u : a(u)— n I 



(3.4) 



then each basis vector x„ £ l 2 (T c i) (v G .Fd) is mapped via n to its abelianization 
Xn £ ^,(Z+) and then II is extended to the whole space t 2 (J-d) via linearity. 
The norm on ^(1/^) is arranged so as to make II a coisometry from ^(Td) onto 
£ 2 with initial space equal to ^^(5^) and with kernel equal to the subspace 

l%{ST d )^ oil 2 (T d ) given by 



(SF^ = <j {/„ W d : fv = for cach n e Z + 

vGFd ■ a.(v)=n 



If we introduce the Z-transformed version II: i/ 2 (JT d ) — > H(kd) via 



n: E/^ u - E 



E /« 

v£jF d : a(t>)— n 



A n , 



then similarly II is a coisometry from T-L 2 {Td) onto Tl(k d ) with initial space equal 
to the subspace 



H 2 {ST d ) := <{ ^ /«(„)«": £ |/„| 2 < oo 



with kernel equal to 



H 2 (sj 7 d) 1 ~ = { f* zV e : E /« = for each n e z + r • 

wG^d u : a(v)=Ti 
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This gives the natural link between the Fock-space norm on formal power series and 
the Arveson-space norm on analytic functions on the unit ball and is the basis for the 
application of noncommutative results to prove commutative results in 6, 251 147) . 

By a commutative (/-dimensional linear system we mean a linear system with 
evolution along the integer lattice Z^ rather than along the free semigroup Td- 
A particular type of such a system is a system of the Fornasini-Marchesini form 
given by 1)1.18(1 . If we specify an initial condition x(0) = a; G X along with an 
input sequence {it°(n)} n6Z d and impose the boundary conditions that x(vi) = 

whenever n is outside the positive orthant Z^ , then the system equations uniquely 
determine a full system trajectory {u(n), x(n), y(n)} consistent with x(0) = x° and 
u(n) = u°(n) for neZf. ' 

If II is the projection map introduced in 1(3.4(1 formally extended to be defined 
on all .^-indexed sequences to generate a Z^-indexed sequence 

II: {u(v)} veJ r d h-> 1 u(v) 

I v : a(v)— n 

H: {x{v)} ve r d ^ < Y x ( v "> 

I v : a(v)— n 

II: {y(v)}ver A ^ I Y y ^ 

I v : a(v)— n 

then one can check the claim: {(IIu)(n), (ILr)(n), (IIy)(n)} ngZ d satisfies the sys- 
tem equations 1(1.18)1 whenever {u(v), x(v), y(v)} V £f d satisfies the system equations 
((1.13(1 . Indeed the first system equation in ((1.13(1 can be rewritten in the form 

d d 

x(v) = Y A k x(k~ l v) + Y B k u(k~ 1 v). 

k=l k=l 

Here we use the convention that 



k v = 



if v = kv' , 
undefined otherwise 



for v a word in Td and k 6 {1, . . . , d} a letter and that x{k v) is interpreted to be 
if fc _1 D is undefined. Summing over v € Td with sl(v) = n then gives 

d d 

(Ux)(n)=Y^k Y x(k~ 1 v) + Y B t Y < k ~^)- 

k—l v : a(v)— n k— 1 v : a(v)— n 

Now observe that 

{k~ 1 v : a.(v) = n} = {v' G T d ■ a(u') = n - e fc } 

and arrive at 

d d 

(ILc)(n) = ^ A fe (n.x)(n - e fe ) + £ B fc (n U )(n - e fc ). 

fe=i fe=i 
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We see that {(IIu)(n), (IIa;)(n), (IIy)(n)} satisfies the first of the system equa- 
tions Hl.lSjl . That {(ITu)(n), (ITar)(n), (IIy)(n)} satisfies the second system equa- 
tion in (|1.18|l is a simple consequence of linearity. Conversely, given a trajec- 
tory {«(n), x(n), y(n)} neZ d of (|1.18(l . let {ui(v)} ve jr d be any ZY-valued .^-indexed 

sequence such that Tlui = u and set a^(0) = x(0). Then the noncommutative 
system equations (I1.13f) recursively uniquely determine a full system trajectory 
{ue,(v), Xi(v), ye(v)} V £F d of (|1.13() with this preassigned input string and initial 
condition. By the claim verified above, it follows that (Hue, Tlxi, Hye) is again a 
system trajectory. By the uniqueness of solution of the initial value problem for 
the system H1.18f) . it follows that {(Hu,Hx,Hy)} = {(u,x,y)}. Thus, any trajec- 
tory {u(n),x(n),y(n)} neIdd can be lifted to a trajectory {ut(v), Xe(v), ye{v)} ve r d 
of (| 1 . 1 H|) - i.e., {ui(v),xi(v),yi(v)} V £r d is a trajectory of (|1.13() such that 

{ILui, Tlx e , Tlyt} = {u, x, y}. 

In this way we view the Fornasini-Marchesini commutative system <|1.18ll as the 
abelianization of the noncommutative Fornasini-Marchesini system 1(1. 13|) . 

Since the commutative Fornasini-Marchesini system (|1.18(l is just the abelian- 
ization of the noncommutative Fornasini-Marchesini system 1(1.13(1 and we have al- 
ready derived the formula 1(2.6(1 for the solution of the noncommutative initial- value 
problem, we see that the solution of the initial-value problem for the commutative 
Fornasini-Marchesini system p. 18(1 is simply the abelianization of the correspond- 
ing formula for the noncommutative case: 

(fig) (A) = C(I - Z(\)A)~ 1 x(Q) + Is (A) • (9u)(A) (3.5) 

where the transfer function Ts(A) for the commutative Fornasini-Marchesini system 
is given by 

T S (A) = D + C(I - Z{\)A)- l Z[X)B. 

This gives a derivation of the transfer function relationship ((3.5(1 (via the connection 
with noncommutative systems) which is an alternative to the usual direct approach 
via commutative multivariable Z-transform (see e.g. |14j1. 
The zero input string simplifies the system to 

( x(n) = Ai^(n-ei)H h A d x(n - e d ) . . 

\ y(n) = Cx{n). 

Given a pair (C, A), we have the option of considering (C, A) as coming from a 
noncommutative or a commutative system. If we consider the associated noncom- 
mutative system, the output string associated with initial state x(9) = x (and zero 
input string) is the 3^-valued function on given by 

O c ,ax = {CA"x} veJ r d 

and (C, A) is considered output stable if this output string is in ly(J-d) for all 
x £ H. We say that the commutative system ((3.6(1 is output stable (and in this case 
we will say that the pair (C, A) is a-output stable) if H(Oc,ax) G l 2 w y{^\.) for all 
x £ H, or equivalently, if HOc,ax) is in the Arveson space Hy(kd) for all choices 
of initial state x £ X. We note that nOc,A^ can be computed explicitly as 

(n<3 c>A z)(A) = C{I - Z{X)A)- 1 x. 
Thus another equivalent formulation of a-output stability is: 
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Definition 3.1. A pair (C, A) is a-output stable means that the function C(I — 
Z(\)A)~ 1 x belongs to Hy(kd) for every x £ Tt, or equivalently (by the closed graph 
theorem), the operator O c A from X to TCy(kd) defined by 

d c>A = n<5 c ,A :x^c(i- z(x)A)- 1 x = J2 [ J2 CAyx A " ( 3 - 7 ) 

is bounded. 

The inverse Z-transform sends the function 

y a (\) = C(I-Z(\)A)- 1 x= ( CA v x\\ n 

neZ^_ \«ea- 1 (n) / 

to the string {y(n)} neZ d with 

y(n) = CA v x, neZf (3.8) 

wGa _1 (n) 

and y 3 - belongs to Hy{kd) if and only if {y(n)} nGZ d £ y(Z^). Thus, the operator 
O c a introduced in Ij3.7|l is the Z-transformed version of the observability operator 

oh, A - J2 CA " x \ ( 3 - 9 ) 



and a pair (C, A) is a-output stable if and only if O c A is bounded as an operator 



from X into £^ y (lA). In this case it makes sense to introduce the observability 



gramian 

Qc,a ■= (O a c, A )*0^ A = (^, A )*% A 
and its representation in terms of strongly converging series 



^A = E 



A*" T C*CA U ) (3.10) 

v,u£a~ 1 (n) 

follows immediately by definitions l|3.9|) , (|3.7|) and the formulas for the inner prod- 
ucts in £l y (1 c l) and Hy(k d ). 

Definition 3.2. We say that the pair (C, A) is a.-observable if Gc a 1S positive- 
definite and exactly ^-observable if Gc a is strictly positive definite. 

By Theorem 12.21 (2) we know that the observability gramian Gc,a satisfies the 
Stein equation l|2.15|l . It turns out that the abelianized observability gramian Gc a 
satisfies a reverse Stein inequality (the reverse of l|2.14|l '). 

Proposition 3.3. Let (C, A) be an el- output- stable pair and let Q c A be the abelian- 
ized observability gramian (|3.10(l . Then Gc a satisfies the reverse Stein inequality 

Qc,a - A\g CjA A x A* d g c>A A d < C*C. (3.11) 

Moreover, the following are equivalent: 
(1) Equality holds in (|3.11l) . 
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(2) A is C-abelian in the sense that 

CA" = CA", whenever v,u € and a(u) = a(it). 



(3.12) 



(3) The observability gramian and the abelianized observability gramian are 
identical: 

Sc,a = Gc.a- 

Proof. It suffices to show that the operator Q given by 

d 



Q:=c*c-g^ A + J2^ A A 3 

3=1 

is positive semidefinite. To this end, plug l|3.1()(l into Ij3.13l) to get 

oo 



(3.13) 



(3.14) 



N=l 



where Qm is given by 
i=i 



E 



m! 



■\& d , : |m|=JV-l 



(N-l) 



— V A*" C*CA V 
1)1 ^ 



E "St E ^ vT C*CA\ 



nGZ^.:|n|=JV u.uea-^n) 

We introduce the notation 

W(n) = J2 A" (n e Z%) 

i(£a _1 (n} 

and extend the notation to the all of Z d by 

W{n) = if n € Z d \Z d + . 
With these definitions we have the equality 

d 

W(n) = 22 W(n - e i )A l for every n = (m, . . . , nd) G Z 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



where e\, . . . , e^ 6 Z^ are defined in i|1.19fl . Write formula l|3.15|l in terms of (|3.16|) 



Q 



N 



£4 

3=1 



E 



mGZ|: |m|=JV-l 



(N-iy. 



W(m)*C*CW(m) 



A, 



£ ^ ■ W(n)*C*CW(n) 

neZ+: |n|=JV 



(3.19) 
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Upon rearranging the terms in the first series in Ij3.19|l and substituting formula 
(|3.18|) into the second series, we arrive at 

= E J2j^§ A J w ( n -^y c * cw ( n - e ^ ( 3 - 2 °) 

neZ^: |n|=JVj'=l 

- E ( E M^( n - ei)*C*CW(n - e,)A, 

{neZ^_: |n|=JV} \*.J=1 

We now consider the terms in i|3.2U[l that correspond to a fixed n = (m, . . . , n<j) S 
(with |n| = N). Denoting the sum of these terms by S n we have 

5 « =E Jffr§ A *i W ( n " ^)*C*CW(n - ej )^ 

" M E ^* M/ ( n - e,)*C*CTU(n - ej )Aj 

= E ^ %")^^(n - e;)*C*CW(n - e,)A,- 

j'=i 

E A*^(n - ei rC7*CM/(n - e,)^, (3.21) 

i,ie{l,...,d} : i/j 

Note that by convention l|3.17[l . the indices i and j in the latter summations vary 
on the set 

Z n = {£e{l,...,d}: n e >0} 
rather than {1, . . . , d}. Furthermore, since 

N -nj = |n| -rtj= n i 

iel„: i^j 

and 

(n - e,-)! = (n - e, - e*)! m (i ^ j), 
one can rewrite the first sum on the right hand side in (|3.21() as 

E n i A *j W (" - e 3 yC*CW(n - e 3 )A, 

= E ^ ~ 6 m n * A *i W( - n - ^rC*CW(n - ej)Aj. 
Plugging this into the right hand side in (|3.21l) leads us to 

^ = £ ^ ~ Cl)l [n*A*W(n - ej yC*CW(n - ej )Aj 
-riiUjA* W{n - ei )*C*CW(n - ej)A 3 ] 
E 2~ M ~ Ix '" - j Ix '" - i ^ 3 - 22 ) 

i,j£X n : i^j 
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where 

R n i j = C [mW(n - ej)Aj - nj W{n - e i )A l ] . (3.23) 

Representation ((3.22|) implies that S n is positive semidefinite and therefore Qn > 
for every N £ N. By 1(3. 14(1 . the operator Q defined in (|3.13|) is positive semidefinite 
which completes the proof of l|3.11|l . 

We now show the equivalence of (1), (2) and (3) in the second part of Proposition 
EP1 

Proof of (1) (2): Assume condition (1), i.e., that the reverse Stein inequality 
113.11( 1 is satisfied with equality. Then representation (|3.22|) implies that R n ,i,j — 
for all neZf. By lj3~2^jl . this means that 

n*CW(n - ej)Aj = n 5 CW{a ~ e 4 )A ( (n £ Z d + ). (3.24) 

Now we shall prove ((3.12(1 by induction (on the length of words v,u £ Td). The 
basis of induction ( \v\ — \u\ — 0) is trivial. Assume that 1(3.12(1 holds true, whenever 
\v\ = \u\ < N. Then in particular, we have for every m £ Zl with |m| < N: 

CW(m)= CA W = ^CA W ° for every w Q £ a 1 (m). (3.25) 

w£a~ 1 (m) 

Now take two words v, u £ of the length N and let 

a(v) = a(u) =: n = (m, . . . , n d ). (3.26) 

If v = v i and u = ui for some v , u £ Td and i £ {1, . . . , c?}, then we have CA V = 
CA U by the induction hypothesis and therefore, 

CA V = CA^'A, = CA^Ai = CA U . 

Let v = vi and u = uj for some i, j £ {1, . . . , d} and i ^ j. By (I3.26|) . sl(v) = n — 
and a(u) = n — ej . By ((3.25(1 , we have 

(N-D) 

CW(n-e 3 ) = ) CA V , (3.27) 

fAT-1)! 

W(n-e.) = ) —CA W . (3.28) 

(n - e t )\ 

Multiplying 13.27(1 and ((3.28(1 on the right by riiAj and njAi respectively, we get 

fAT-1)! - (N - 1)\ ~. (N - 11! 

mCW(n - e,)A 7 = m) ttCAM,- = (-CA^ = mnA r -LCA v 

(n-e,)! ' (n-e,)! n! 

and 

(AT-1)! - (AT-ll! ~. fAT-ll! 

(n-ei)) J (n~ ei y. n! 

By ((3.2411 , the left hand side expressions in the two latter equalities are equal. Upon 
comparing the right hand side expressions we get CA V = CA" , i.e., A is C-abelian 
as wanted. 

Proof of (2) ==> (3): Assume now that A is C-abelian, i.e., that ((3.12|) holds. 
Then the identify Q^, A = Gc,a is an immediate consequence of the series represen- 
tations 1(3.10(1 and 1(2.10(1 for Gc.a an d Gc,a respectively. 

Proof of (3) => (1): We know from Theorem EH (2) that Gc,a satisfies the 
Stein equation, i.e., Gc,A satisfies the Stein inequality 1(3.11(1 with equality. Hence 



rijCW{n - ei)Ai = ^CAM* = n \- '-CA ul = nj rn - — — — — YA U . 
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trivially Gq a satisfies 1)3. with equality whenever Gq a = (7c, a- This completes 
the proof of Proposition EH D 

Example 3.4. If (C, A) is an output-stable pair, then by Theorem 12.21 (2) Gc,a 
satisfies the Stein equation (|2.15|l and hence in particular 

Gc.a - AIGc.aAi A* d GcMd > 0, 

We now show that, for the abelianized case, the inequality in the reverse Stein 
inequality satisfied by the abelianized observability gramian Gq a can ^ e strict hi 
the strong sense that the quantity Gq a~~ A^Gc aAi — ■ ■ ■ — A^Gq aAcL is not even 
positive semidefinite. As an example, let 



C = [1 0] , 

A straightforward calculation shows that 

5c',A-AlG^ A Ai-...-A* d G^ A A 

which is not positive semidefinite. 

Condition 1)3.1 2JI is worth a formal definition. 

Definition 3.5. Let C £ C{X,y). A d-tuple A = {A u . 
tors on X will be called C-abelian if 13.12fl holds. 





" 


1 

2 


0" 




"0 





1" 

2 


At = 











, A 2 = 


1 

2 










1 

. 2 




















, Ad) of bounded opera- 



One obvious way for a given operator d-tuple A to be C-abelian is for A itself 
to be commutative, i.e., for AiAj = AjAi for all 1 < i,j < d. We next show that, 
under an observability assumption, this is the only way. 

Proposition 3.6. Suppose that the output-stable pair (C, A) is observable and that 
A is C-abelian. Then the d-tuple A is commutative. 

Proof. Since A is C-abelian, relations (|3.12|l hold. Fix i,j e {1, . . . , d} and note 
that by (EHJ), 

CA v AiAj = CA mj = CA VJl = CA v A,jA % for every v e T d , 
since a.(vij) = &(vji). Thus, 



CA v (A i A j 



AjAi)x 







for every v £ Td and x S X. Since the pair (C, A) is observable, we have by (|2.13(l 

(^A-^A-j — Aj A-i)x = 
holding for every x £ X, which proves the commutativity relations 



A.;Aj 



AjAi 



for i,j = 1, 



□ 



and completes the proof. 

Corollary 3.7. Suppose that (C, A) is an observable output-stable pair. Then the 
abelianized observability gramian coincides with the observability gramian 

Gc,a = £c,a 

if and only if the operator d-tuple A is commutative. 
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Proof. Combine (2) (3) in Proposition 13 . HI with Proposition 13. 61 



□ 



We next show that the observability gramian always dominates the abelianized 
observability gramian. 

Proposition 3.8. Let (C, A) be an output-stable pair. Then: 

(1) (C, A) is also a- output- stable with 

Qc,a < Sc.A- (3.29) 

(2) Equality occurs in (|3.29|) if and only if A. is C-abelian: 

CA V = CA U whenever u,v £ Td with a(u) = el(v). 

Proof. Note that the second statement in Proposition 13.81 is just a restatement of 
(2) (3) in Proposition 13. 31 Thus it suffices only to prove the first statement. 

By definition, output-stability of (C, A) simply means that Gc,a is bounded, 
while a-output stability means that Qq a is bounded. The fact that a-output sta- 
bility follows from output-stability therefore follows immediately from the general 
inequality (|3.29l) . Thus it suffices to prove l|3.29|l . For this purpose, recall that 

(Qc,ax,x) = \\O c ,Axf ey{:Fd) = £ \\CA v xf y , 

while 



{Gb,A.x,x) = \\Ob,A* 



^ Inl 



E CA% 

Jj£a -1 (n) 



By the Cauchy-Schwarz inequality we have 



•u£a -1 (n) 

Therefore 




CA v x\\y 



* E 

?)£a _ 1 ( n ) 



CA V x\\y 



n ! 



CO; 



C,A J 



< 



y n! 

E E 



CA 11 

^Ga -1 (n) 

\CA* X \\y 



x\\y 



\0 C ,AX\ 



and l|3.29|l follows as wanted. 



□ 



Example 3.9. The converse of Proposition I3T81 part (1) can fail, i.e., there exists 
an output pair (C, A) which is a-output- stable but not output- stable. For example 
take 

C = [1 0] , At 

Then C(I - AiAi - = [l 2Ai 2A 2 ] . Hence 

Oc,a -x^C{I- A1A1 - \ 2 A 2 y 1 x 



' 


2 


0" 




"0 





2" 











, A 2 = 


1 








-1 
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maps X — C 3 into %{k%) and thus (C, A) is a-output stable. To show that (C, A) 
is not output stable, note that 

"2" 

(AiA 2 ) n =00 
(-2) 1 

and therefore, C{A x A 2 ) n = [ 2 n ] , so that for x = [ 1 ] T , 

£ ||CA^>^2»=oo 

v£j-d n>0 

and therefore, the pair (C, A) is not output-stable. We conclude that a-output- 
stability has no obvious characterization in terms of positive semidefiniteness of 
some solution of a Stein inequality as in the noncommutative case (see Theorem 

E2(2)). 

As a corollary of the gramian inequality (|3.29|l in Proposition 13.81 we have the 
following. 

Corollary 3.10. Let (C, A) be an output-stable pair. Then: 

(1) KerC^A C KeiGcA- Hence, if (C, A) is ^.-observable (respectively, ex- 
actly a.- observable, then (C, A) is also observable (respectively, exactly ob- 
servable). 

(2) The subspace Ker Gc, a = Ker Oc, a is invariant under the operator Aj for 
each j = 1, . . . , d. 

(3) The subspace Ker Gc A is invariant under Aj for each j — 1, ...,d if and 
only if Ker Gc,a = Ker Gc, a • 

Proof. Statement (1) is an immediate consequence of the inequality Ij3.29|l . State- 
ment (2) is easily checked from the definition of Oc,a- Sufficiency in statement 
(3) is then a consequence of statement (2). It remains only to verify necessity in 
statement (3). 

Assume therefore that KerC/J A is invariant under Aj for each j = 1, . . . , d. Let 
x be a vector in KerQ^ A . Then by the assumed invariance, A u x £ KerQ* A for 
every u £ T&. Then we have 

Ker CA vu x = for every u £ T d . 

ti£a" 1 (n) 

Then letting n = we get CA u x — for every u £ Td and therefore, x £ Ker Gc,A- 
Thus, Ker Gq a ^ Ker Gc,a and since the reverse inclusion holds by the first state- 
ment, equality follows. □ 

Example 3.11. We observed in part (1) of Corollary 13 . 1 01 that a-observability for 
an output-stable pair (C, A) implies observability. We now give an example to show 
that the converse can fail, i.e., there exists an output-stable observable pair which 
is not a-observable. For this purpose, let d = 2, X = C 4 , y = C, C = [0 l] 
and A = (Ai,A 2 ), where 
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Then the pair (C, A) is output stable. Now we show that (C, A) is observable but 
not a-observable. Indeed, since 



we have 



and 



CA X 



CA1A2 — [ lis 2S6 2B6 °] 



1[0 -1 1], CA 2 = -±[1 1] 



CA 1 A 2 



256 L 



1 0] 



Now it is clear that KerC n KevCAi n KerCL4 2 n KerCAiA 2 = which implies 
that f] ve jr d Ker CA V ~ 0. Therefore, the pair (C, A) is observable. To show that 
(C, A) is not a-observable we first compute 



/ - A1A1 - X 2 A 2 



' 16 

Ai _i_ Aa 

16 \ 16 
A2. 

A 16 

A2 

16 



_Al 

\ 16 
1 _ Ai _|_ 
1 16 ^ 
A2 
16 





16 








Ai 


+ Aa 


16 




1 


A 1 / 






A, 16 




16 



Aa 

_ Ai _i_ 

A2 
\16 

1 _ Ai _ 

1 16 



Aa 

16 



A2 

16. 



A straightforward calculation gives 

<2(Ai,A 2 ) := det (J - - 
_ Xi XiX 2 



X 2 A 2 ) 
XfX 2 



\2 
A 2 



A 1 A 2 



A X A^ 



\4 
A 2 



16 128 2048 64 2048 16384 16384 
Note that 

[2/1 2/2 2/3 2/4] := C(7 - X1A1 

is the bottom row of the matrix (I — X\Ai + X 2 A 2 
adjoint formula for the inverse of a matrix to get 



A 2 ^ 2 r i 

1 and we use the standard 



272 



1 



d(Ai,A 2 ) 

Then it follows that the nonzero vector x 



_Aa 

\ 16 
A 2 

16 



Aa 

16 



_Al 
Aa 6 

16 







1 A2, 

Ai 16 

16 



[0 1 0] T satisfies 



C(/-AiA 1 -A 2 A 2 )- 1 x = 
and therefore, the pair (C, A) is not a-observable. 

3.2. Observability-operator range spaces and reproducing kernel Hilbert 
spaces: the commutative-variable case. We seek the analogue of Theorem 
ll.2l for the commuting multivariable case. We extend multivariable power notation 
(|3.1(l to any <i-tuple A = (Ai, . . . , Ad) of commuting operators on a space X: 

A n := A^A^ 2 ■ ■ ■ A n d d . (3.30) 

Note the connection between the commutative powers A n (with n £ Z^) and the 
noncommutative powers A" (with v € J-'d) in case A is a commutative operator 
d-tuple: 

N\ 



A" = A n where n = a(v), ^ A v XA V = ^A.'"XA' 



vGF d : \v\=N 



nSZ° : |n|=JV 
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for any operator X on X . In case (C, A) is an output stable pair with A a commu- 
tative operator d-tuple, the formulas (|3.7[) . I|2.1()|l and (^.lOjl for Oq A , Gc,a and 
Gq a collapse (in view of (|3.2|0 to 



Oc,A- x^C(I - Z(X)A)- 1 x= ^y(CA n a;)A n (3.31) 

and 

Gc.a = Gca= Y ^A* n C*CA n . (3.32) 
c — ' n! 

nGZj 

We next observe that a natural commutative counterpart of operators Sj intro- 
duced in 12.25fl are the operators M Xj of multiplication by the coordinate functions 
of C d for j = 1, . . . , d acting as contractions on the Arveson space Tiy(kd). We will 
call the commuting d-tuplc Ma := (M\ 1 , ■ ■ ■ ,M\ d ) the shift of Hy(kd), whereas 
the commuting rf-tuple M A := (M Xi , • ■ ■ , M\ d ) consisting of the adjoints of M\ j 's 
(in the metric of Tty(kd)) will be referred to as to the backward shift. Recall that 
monomials A n form an orthogonal basis for TC(kd). As we have seen, 

<A n , X m )n( kd) = { H! if n = m (3.33) 
I otherwise. 

A simple calculation based on (|3.33(1 gives 



id 



M* x .X m = ^X" 1 -^ {m 3 > 1) and M x .X m = (rrij = 0) (3.34) 



m 



where m = (mi, . . . , m^) and ej 6 U\_ is defined in Q1.19|l . More generally, 



m!|m-n|! „ 

■ A ,ii wij > Tij tor j = 1 , . 



(M A ) n A m = ^ |m|!(m-n)! ' 3 ~ 3 J (3.35) 

0, otherwise, 

where according to l|3.30|l 

(Mir ■■= (M* Xl r (mix 2 ■ ■ ■ (Mir ■ 

The following proposition includes the analogue of Proposition 12 . 91 for the present 
commutative setting. 

Proposition 3.12. Let M A be the d-tuple of backward shifts on TLy{kd) and let 
G: TLyikd) — > y be the operator of evaluation at 6 M d 

G: /(A)-»/(Q). (3.36) 

Then: 

(1) For every f e "Hy(kd) and every X € M d we have 

d 

f(X)-f(0) = J2\ j (M* x .f)(X). (3.37) 

3=1 

(2) The pair (G, M A ) is isometric: 

I - M Xl M* Xi - ... - M Xd M* Xd = G*G. (3.38) 
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(3) The abelianized observability operator associated with the pair (G, M^) is 
the identity operator: 



(4) The d-tuple is strongly stable, that is, 



(3.39) 



Nl, 



uS^d : |«|=iV 

for every f GHy(k d ) 



nGZ^: \n\=N 



(3.40) 



Proof of (1): One can easily verify the identity (|3.37(l on monomials y ■ A m (with 
y € y and m € Z^_) using IpOIj) . Then the result follows for all / 6 U y {k d ) by 
linearity and continuity. 

Proof of (2): Note that G* : y — > TLy{k d ) is the identification of a vector y € y with 
the constant function y € TLy{k d )- We then see that H3.38|l is simply the operator 
expression of l|3.37(l . 

Proof of (3): From lEHSjl and (|3~^|) we see that 

G(M* x ) n f=^-f n if /(A) = ]T / m A m and n e Z% 

and therefore, according to definition l|3.31|l . 

GM1 f:= E ^(G(Ml) n /) A n = ^ / n A n = /(A). 

Since the latter equality holds for every / e l~ty(kd), 13.39|l follows as asserted. 

Proof of (4): This can be derived directly from (|3.35|) or via Proposition 12 . 71 since 
Oq m , = I and therefore, the pair (G, M^) is exactly observable. □ 

Remark 3.13. Note that in contrast to the noncommutative case (Proposition 
the operator 



R = 



ML 



M 



G 



Hy{k d ) (Hy(k d )) d ® y 



~P 








"(^(fe<i)) d " 







y 




y 



is not unitary (just isometric). A simple calculation shows that 
/ - RR* = 

u uj L y 

where P is the orthogonal projection of {'Hy(k d )) d onto the subspace 



e (^(fc d )) d : 5>A W = o 





V 
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If a pair (C, A) is a-output stable, then the observability operator Oq a : X — * 
7~Ly(kd) is bounded and its range 

RanOg, A := {C(7 - Z(\)A)~ 1 x : x <E X} (3.41) 

is a linear manifold in TCy(kd). We have the following partial analogues of part (3) 
of Theorem PI 

Theorem 3.14. Let (C, A) be an a-output stable pair. Then: 

(1) RanOJ A with the lifted norm 

\\C(I - ^(A)^)- 1 ^^ a) = \\Q*x\\ x (3.42) 

where Q a is the orthogonal projection of X onto (Ker^g. A )- L is isometri- 
cally equal to the reproducing kernel Hilbert space TL{Kq A ) with reproducing 
kernel Kq A (A, £) given by 

#c,a(A, = C(J - ^(A)A)- X (/ - ^Z(C)*)- X C* (A, C G M d ). 

(2) Ran Op A with norm inherited form TLy{kd) is a reproducing kernel Hilbert 
space 7~i(K c ' A ) with reproducing kernel K c ' A (A, £) given by 

k*--(\,c) = c(i - z(X)A)- 1 (g^ A )- 1 (i-A*z(cyr 1 c* (A,C G B d ). 

We next discuss separately the case where A is C-abelian and then the general 
case. 

3.2.1. T~C {Kq a ) for the case where A is C-abelian. In case (C, A) is an a-output- 
stable pair with A C-abelian, then we have the following commutative analogue of 
Theorem |2~TU1 

Theorem 3.15. Let (C, A) be a contractive a.- output- stable pair such that operator 
d-tuple A is C-abelian. Then: 

(1) The intertwining relations 

MZ j d% tA = dZ iA A j forj = l,...,d (3.43) 

hold, and hence the linear submanifold Ran Ofj A ofTLyikd) is Mt-muan- 
ant. 

(2) The operator Oq a maps X contractively into TLy{kd). This mapping is 
isometric if and only if (C, A) is isometric and A is strongly stable. 

(3) If M := RanOg. A is given the lifted norm l|3.42[) (so M is isometrically 
equal to H.(Kq a ) by Theorem \S.14\ (1)), then the difference- quotient in- 
equality 

Ell^/llW? A )^ll/llWe A )-|l^ )l^ 

i=i 

holds for every f G Tl(K^, A ). Moreover, the difference- quotient identity 
EP^/llW ?A ) = ll/llWe A )-|l- f (°)l^ 
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holds for every f 6 Ti.(K^ A ) if and only if the subspace (Ker<?c,A) * s 
A-invariant and the restriction (C , A ) (defined in 12.53J1 ) 0/ (C, A) to 
i/ie subspace (Kei Qc. a) 1 ' is isometric. 

Proof. By (|3.31[) and (|3.34|) . we have for every 2 e A", 
(Af Xj )*5& iAa: = (M A .)* 




and l|3.43|l follows. This completes the proof of statement (1) in the theorem. 
Since the pair (C, A) is contractive and A is C-abelian, we have 

Gc,a = Gc,a <Q = Q a <L 

Therefore, 

\\0^ A x\\ ny (k d ) = (eh, A x, x)\ < \\Q*x\\ x = \\d^ A x\\ n{KbA) < \\x\\ x . (3.44) 

Now the arguments used in the proof of Theorem 12.101 can be used to prove the 
remaining statements in the theorem. □ 

For the converse direction we have the following result. 

Theorem 3.16. Let Ai be a Hilbert space of y -valued functions included into 
TLy{kd) and let us assume that Ai is M^- invariant. 

(1) If the inequality 

d 

EH M AVIlA4<ll/llA4-|l/(0)g (3-45) 

i=i 

holds for every f £ Ai, then Ai = Ran(0ji, A ) for a contractive and exactly 
observable (with respect to Ai) pair (C, A) with the commutative d-tuple 
A = (Ai, . . . ,Ad). In particular, Ai is contractively included in T~Ly(kd). 

(2) If the equality 

d 

EH m a 3 /IIa4 = II/IIa4- ll/(°)Hy ( 3 - 46 ) 

holds for every f £ Ai, then Ai = Ran(OJ A ) for an isometric and exactly 
observable (with respect to Ai) pair (C, A) with the commutative d-tuple 
A. By part (1), Ai is contractively included in Hy(kd)- Moreover, it 
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is isometrically included in Hy{kd) if and only if the restriction of the 
backward shift to Ai is strongly stable, i.e., 

lim £ ^j\\(M* x rf\\ 2 M =0 for every feM. (3.47) 

nGZ^ : \n\=N 

Proof. Define operators C : M —> y and Ai,...,A d : M — » M by 

C = G|a4: /^/(Q) and A 5 = M* X .\ M (j = l,...,d). (3.48) 

Thus, the (i-tuple A is the restriction of the backward-shift tuple to Ai. By 
part (3) of Proposition ^. 121 it follows that <D^ A = Im an d thus, the pair (C, A) is 
exactly observable (with respect to Ai) and the range of the associated observability 
operator Oq a coincides (algebraically) with Ai . Now we write (|3.45j) in terms of 
operators l|3.48|l as 



^\AifTM<\\ftM-\\ c fty (/ e -M) 
.7=1 

and conclude that the pair (C, A) is contractive. Similarly, assumption 
means that the chosen pair (C, A) is isometric. Furthermore, if Ai is included in 
7~Ly(kd) isometrically, relation (|3.47|) holds since is strongly stable (see part 
(4) of Proposition 13 . 1 2)l . Conversely, if (|3.47|l holds, that is, if the commutative 
rf-tuplc A = (Ai, . . . , Ad) is strongly stable on Ai, the Stein equation (|2.15|l has a 
unique positive semidefmite solution. Since the pair (C, A) is isometric (recall that 
we are proving isometrical inclusion under assumption H3.4fifl ). this unique solution 
is the identity operator. On the other hand the observability gramian Qq a = Gc,A 
defined by the convergent series (|3.32() satisfies the same Stein equation (as observed 
in part (2) of Theorem 12.2(1 . Thus, Gc a = Sc,a = I- Note that the inequality 
(l2~T%j) holds with H = g c ,A = I, i.e., 

A *v T C *CA v <I- A* yT A v . 

ver d : \v\<N ver d : \v\=N+l 

Taking strong limits as N — > oo and noting that / — Qc,a — J2veF d C*CA V 
then gives 

I < I — s-limA^oo Y A* yT A v 
ver d -. \v\=n 

from which the strong-stability of A follows. Then Ai — Ran(Op A ) is isometrically 
included in Tty(kd) by statement (2) in Theorem 13. 151 □ 

We have the following analogue of Theorem 12.131 for the present commutative 
situation. 

Theorem 3.17. Suppose that (C, A) and (C,A) are two observable output-stable 
pairs with both A and A commutative such that A {\, £) = if- ^(A, £) for all 

A, £ 6 M d . Then there is a unitary operator U : X — > X such that 

C = CU and Aj = U^AjU for j = 1, . . . , d. (3.49) 



44 



J. A. BALL, V. BOLOTNIKOV, AND Q. FANG 



Proof. Suppose that (C, A) and (C, A) are as in the hypothesis of the theorem. The 
identity of the kernels Kq a and K~, ^ implies equality of the respective coefficients 

of A n C m for each n, m e Z f {: 

|n|! |m|! Inl! Iml! — ~ ~ 

CA n A* m C* = I CA n A* m C*. 



n! m! n! m! 

If we define a mapping U by 

U: A* n C*y h-> A* n C*y (3.50) 
it follows that U extends by linearity to an isometry from 

V := span{A* m C** 2 /: m e %% and y e y} 

onto 

K := span{A* m C7*y: m e 1% and y e y} 

Since both (C, A) and (C, A) are observable, we see that T> is dense in X and 
that 1Z is dense in X. Hence U extends to a unitary operator from X onto X by 
continuity. From the defining equations (|3.5U[1 for U we see that 

UC*=C* and UA*=A*U. 

By taking adjoints and using that U is unitary, we arrive at the intertwining equa- 
tions l|3.49fl as wanted. □ 

Theorem l3 . 1 71 can be adapted to give the following result concerning containment 
between two backward-shift-invariant subspaces rather than equality; the finite- 
dimensional case appears as Proposition 1.2 in |19| . 

Theorem 3.18. Let M. and M. be two backward- shift- invariant subspaces of the 
Arveson space Tiyikd) with realizations 

M = R&nd c ,A and _M=Ran(5 5A (3.51) 

where the d-tuples A — [A\, . . . , Ad) £ X d and A — (A\, . . . , Ad) £ X d are com- 
mutative and strongly stable and the pairs (C, A) and (C, A) are isometric. Then 
A4 C M. if and only if there exists an isometry V : X — > X such that 

C = VC and VAj = AjV (j = l,..., d). (3.52) 

Proof. The necessity part is clear. For the sufficiency part, assume that M. C ftA. 
By Theorem 13. 151 there exist unitary operators U : M — * X and U : M. — > X such 
that 

U*A,U - M* x . \ M , U*A,U = M* x . \ M (j = 1, . . . , d). 

and 

CU = G\ M , cu = g\^ 

where the operator G : 7iy — > y is defined in (|3.36|l . Let T: M. — > M. be the 
inclusion operator. Clearly T is isometric. Then the operator V = U*IU : X — > X 
is isometric and satisfies (13.521). □ 
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3.2.2. H(Kq a ): The general case. In case the a-output-stable pair (C, A) is such 
that A is not C-abelian, it can happen that the associated reproducing kernel 
Hilbert space is not invariant under the backward-shift tuple M^, as the following 
example shows. 



Example 3.19. Let 

C = 



yj o ' 

2 u 

o 4 





h o 



\ 




Then a straightforward calculation gives 

^c,a(A,C) - C(I-Z(X)A)- 1 (I-A*Z(Crr 1 C* 



(4-AiA 2 )(4-CiC 2 ) 



'2 A 2 " 




" 2 Ci " 


Ai 2 




. C 2 2 _ 



Thus A (A, iu) is positive definite on " 
by the two rational functions 

4 



and the space H(K^, A ) is spanned 



/i(A) 



4- AiA 2 



2 
Ai 



and / 2 (A) 



4 - AiA 2 



A 2 
2 



Furthermore, since 



and since 



it holds that 



4Ai 



ni 



ni + n 2 



m — 1 \n 2 
2 > 



4 - AiA 2 



00 P+ 1 P 



4J 



4Ai 



4- AxA 2 



E 

J=0 



i + 1 /AiA 



2.7 + 1 



The latter function is rational if and only if the single- variable function F(z) = 
Y^jLo 2j+i zJ ^ s ra tional. By the well-known Kronecker theorem, F in turn is ratio- 
nal if and only if the associated infinite Hankcl matrix 



= s 



[si+j] t ~ =0 where s k = 



k + l 



2k + 1 

has finite rank. However one can check that the finite Hankel matrices H„ = 
i s i+jYij=o nave f uu rank for all n = 0,1,2,... and hence F(z) is not rational. 
Therefore, M% fi does not belong to H(Kq a ) and hence H(K* A ) is not invariant 
under M Xi . 

For the general case, there is a simple replacement for M.* X \ H ( K * A y Specifically, 
given an a-output-stable pair (C, A), wc define an operator-tuple T = (T 1; . . . , T d ) 
on Ran Og, A by 



TjO, 



Oc A AjX for x e (K.erGc,A) ± and j = 1, ■ ■ ■ ,d. 



(3.53) 
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Wc then have 

/(A)-/(0) = C(I - Z(X)A)- 1 x - Cx 

= C(I - Z{X)A)- 1 Z{X)Ax 
d 

= Yj X i C ^ ~ Z{\)Ay 1 A :jX 
i=i 
d 

= ]TV( a c,AA^)(A) 

i=i 

d d 

= ^A,-(T j a^)(A) = ^A,-(T i /)(A). (3.54) 
i=i i=i 

We next give the following analogue of Theorem 13 . 1 51 for the general case. 

Theorem 3.20. Let (C, A) be a contractive pair with C £ £(X,y) and A = 
(Ax,...,A d ) £ C{X) d . Then: 

(1) The Z -transformed observability operator Oq A is a contraction of X into 
the reproducing kernel Hilbert space TL{Kq a ). It is an isometry if and only 
if the the pair (C, A) is a- observable. 

(2) The space Tt(K^ A ) is contractively included in the Arveson space r Ky{kd); 

it is isometrically included in TLy{kd) if and only if A (as an operator 
from X into TLy{kd)) is a partial isometry. 

(3) For every function f £ 7i{K^, A ) it holds that 

d 

/(A)-/(0) = ^A i (T i /)(A) (AeB d ) (3.55) 
3=1 

and 

d 

E II^/H«(^ A ) < II/IIWsa) - ||/(0)|1 ^ (3 - 56) 

3=1 

where T%, . . . ,T<j £ C(H(Kq a )) are t/ie operators defined in (|3.53|l . 

(4) Equality holds in (|3.56[) /or every / G H(Kq a ) i/ and on/?/ i/ i/ie subspace 
(Ker ^ca) 1 is A-invariant and the restriction (C°, A ) (defined in (|2.53|) ) 
o/(C, A) to £/ie subspace (KerC/cA) * s isometric. 

(5) J/ H{Kq a ) is isometrically included in TLy{kd), then Tj = M£. for j — 
1, . . . ,d and therefore, TCy(kd) is 'M.* x -invariant. 

Proof. Since the pair (C, A) is contractive, the identity operator H = Ix solves the 
Stein inequality (E33J). Then ££ A < £fc iA < i# (by part (1) of Proposition El 
and part (2) of Theorem |2~21) . Thus, 

<5S,a < Q a < J* 

where Q a is the orthogonal projection of X onto (Ker Qq A ) . Therefore it holds 
for every x £ A" that 

II^,a^II^(M = (fic^. *>l < IIQ a *ll* = I|0&,az|I«(* S<a ) < IMU- (3-57) 
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We have the equality instead of the first inequality in (|3.57|) if and only if A = Q a : 

that is, if and only if Oq A is a partial isometry. Furthermore, the second inequality 
in (|3.57|) can be replaced by equality if and only if Q & = Ix, i.e., if and only if the 
pair (C, A) is a-observable. This completes the proof of the two first assertions in 
the theorem. The multivariable difference-quotient relation (|3.55|1 follows by the 
calculation l|3.54fl . Furthermore, for every x G (Ker^J a) ' 

ll<5c,A*lk* &;A ) = IIQ a *ll* = MW 

and thus, Oq a maps unitarily (KerQ^j^) 1 - onto H(Kq A ). Therefore, by l|3.53|l . 
Tj is unitarily equivalent to the compression of Aj to (Ker Q^, A ) ± and hence 

Pi|| < ||A,-|| forj = l,...,d. 

In particular, Tj G C{TL{K^ A )). For an element / — A x G H(K^ A ), we have 

d d 

E \\Tjf\\n(K£ A ) = X H^i^&.A^IIwCKg A ) 
3=1 3=1 

= X] II^C , ,A-4' af llw(K5 A ) 

3=1 
3=1 

d 
3=1 

<H^-ll^ = ||/||^ (KS|A) -|]/(o)g 

where the first inequality holds since Q a < I and the second since (C, A) is a 
contractive pair. This proves inequality (|3.56|l and it is readily seen that equalities 
hold throughout in the last calculation for every x G (KerC/g. A )- L if and only the 
subspace (Ker Qq a )~*~ is A-invariant and the restriction (C°, A ) (defined in (|2.53|l ) 
of (C, A) to the subspace (KerCJg, A )- L is isometric. 

Finally suppose that H(K* A ) is included isometrically in Tiy(k ( i)- Then the 
assumption l|3.56fl becomes 

d 

E ll T i/ll« y (*«.) < \\f\\H y (k d) - \\f{0)\\y f or every / G M. (3.58) 

3=1 

Then we take the inner product of both parts in equality l|3.54|l with /: 
(/-/(0),/)« y(fcd ) = 11/11^) -||/(o)||5, 

and 

d d 

^2(M Xj Tjf, f) H y(k d ) = ^2( T jfl M * x .f) Hy(kd) . 
3 = 1 3 = 1 

Thus, 

d 

\\f\\Hy{k d ) - \\fm\y = ^(Tjf, MlJ) Hy(kd) . (3.59) 

3 = 1 
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For any / in Tiy{kd), applying the identity 1)3. 38|) to / and then taking the inner 
product with / gives us 

d d 

\\f\\H y {k d ) - \\fm\l = £<M Ai My, f) ny[kd) = \m 3 f\\ 2 nyikd) . (3.60) 

Now we conclude from (|XHU|) . j53Hl and [I33gjl that 

d 

ll/lltt^*,,) " ll/(0)|& = EH^/UwyOW 

3=1 

d d 
3=1 3 = 1 

from which we get 

d d 

3=1 3=1 

d d 

o > -E( T ^' M A 3 />« y (M+Eii r ^ii^(M- 

3 = 1 3 = 1 

d 

Adding these inequalities and using that 'S^(Tjf, M^.f) is real then gives 

3 = 1 

d 

> Yl (ll M A,/f - (T,f, M* x .f) - (M* Xj f, Tjf) + \\Tjf\ 

3 = 1 
d 

= Eh m a j /- t 3/IItW 

3 = 1 

Therefore, M x .f — Tjf for j = 1, and for every / e yV( as asserted. This 

completes the proof of Theorem 13.201 □ 

3.3. The Gleason problem: a uniqueness result. Let Ad be a Hilbert space 
of y~ valued functions. A tuple T = (Ti, . . . . Tj) of operators Tj G Ad is called 
a solution of the Gleason problem (see j^OH]) if relation l|3.55fl holds for every 
f E Ad. Let us say that T is a contractive solution of the Gleason problem if in 
addition 

d 

£\\T,f\\M<\\f\\%l-\\f(P)\\y every f E Ad (3.61) 

3 = 1 

or, equivalently, if the pair (T, G) is contractive where G : Ad — >• y is defined by 

G: /(A) -►/(()). (3.62) 

We have the following analogue of Theorem 13.161 characterizing contractively in- 
cluded subspaces Ad of Hy(kd) of the form Ad = TL(K^, A )\ for the general case 
where Ad is not M^-invariant, one simply replaces with some contractive so- 
lution T of the Gleason problem on Ad. 
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Theorem 3.21. Let Ai be a Hilbert space of y -valued functions and let us assume 
that there exists a contractive solution T = (Ti,...,Td) of the Gleason problem 
(i.e., Tj S C(M) such that and (f5~STfl hold for every f G Ai). Then Ai is 

isometrically equal to a reproducing kernel Hilbert space Ti.(K^, A ) for a contractive 
pair (C, A). Therefore, Ai is contractively included in the Arveson space TLy{kd). 

Proof. Take C = G\ M where G is given by (H02J), A = T on Ai. Then IpHTT) says 
that (C, A) is contractive. Iteration of l|3.55|) says that, for each / e Ai, 



/(A) = X h 



ii=i 



(T n f)(0)+J2 X n 



(T h T h f)(p) + ]T X h [(T j3 T j2 T h f(Q) 



J3=l 



■■■+E ^^•••r i2 T J1 /)(o) + •••]•••] 

3fc=l 

This unravels to the tautology 

f(X) = C(I-Z(X)A)- 1 f 

so we recover Ai as A-l = Ran0g, A with — Z(')A) _1 /llw(A£ A ) = i.e., 
Ai = H(Kq a ) isometrically. From the fact that (C, A) is contractive, we have 
seen that Gc a — Sc,a < ^A4- Then 

||C(/-Z(.)A)- 1 /||^ (fed) - {Qh,Af,f)M < \\f\\ 2 M 

and we also have the contractive inclusion property. □ 



Combining Theorems I3.2UI and 13.211 gives the following uniqueness result for 
contractive solutions of the Gleason problem on a subspace M. contained in TLy{kd) 
isometrically. 

Theorem 3.22. Suppose that M. is a subspace of y -valued functions contained 
in TLy{kd) isometrically and that T = (Ti, . . . ,Td) is a contractive solution of the 
Gleason problem on Ai. Then Ai is M^- invariant and T = M^. 

Proof. By Theorem l3.21l there is a contractive pair (C, A) so that Ai = TC(K^ A ) 
isometrically. As Ai is contained in Hy (kd) isometrically, we conclude that Ti(K^ A ) 
is contained in TLy(kd) isometrically. Part (5) in Theorem l3 . 2Ul then asserts that the 



subspace Ai — H(Kq A ) is M^-invariant and that Tj 



Ml forj = l, 



□ 



We note that the proof of Theorem 12.131 is like the proof of the State-Space- 
Isomorphism Theorem for structured noncommutativc multidimensional linear sys- 
tems in |12j . It is known that the State-Space-Isomorphism Theorem (and related 
Kalman reduction procedure) fails in general for commutative multidimensional 
linear systems — see e.g. [201 f° r a recent account of the situation. The fact that 
uniqueness does hold in the special commutative situation in Theorem 13.171 shows 
that the technique in the proof of the State-Space-Isomorphism Theorem is sal- 
vageable in special commutative situations. 

A uniqueness result for solutions of the Gleason problem somewhat different 
from that in Theorem 13.221 was obtained in [3]; rather than assuming that T is a 
contractive solution of the Gleason problem on Ai — H{Kc,a) contained isometri- 
cally in TCy(kd) as in Theorem l3.22l Alpay and Dubi in 3 assume instead that T is 
a commutative solution of the Gleason problem and are then able to conclude that 
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necessarily T = M^|». This latter result can be seen as an immediate consequence 
of our Theorem 13 . 1 71 above since, by the construction in the proof of Theorem 13. 2 II 
solutions (C, A) of A = K are in one-to-one correspondence with solutions T 
of the Gleason problem. We illustrate the preceding analysis by two examples. 



Example 3.23. Consider the subspace Ai 
the operators T 0j i and T Qj2 on M by 

T a ,i ■ f i-> P + aa\ 2 , 



span{l, Ai, A2} C 7i{k2) and define 
T a ,2 ■ f i->- 7 - aa ^i (3.63) 



where /(A) = a + /3Ai + 7A2 is the generic element in Ai and where a is a fixed 
complex number. It is readily checked that 

/(A) - /(0) = /JAi + 7A2 = Ai(T Qjl /)(A) + A 2 (T Q , 2 /)(A), 

so the tuple (T a ^,T a>2 ) solves the Gleason problem on Ai. Let A a ,i and A a ^ be 
the matrices of T a .i and T a2 with respect to the basis {1, Ai, A 2 } of Ai and let C 
be the matrix of the operator G : Ai — > y defined in H3.62(l : 

"0 1 0] I" 1" 

C = [1 0] , A a . x = 0, A a , 2 = -a . (3.64) 

[a oj [ 

A straightforward calculation shows that 

C(I - \!A aA - \ 2 A a>2 )- 1 = [1 Ar A 2 ] 

which realizes Ai as the range of the observability operator of a pair (C, A a ). 
Different choices of a in (|3.64(l lead to non-equivalent realizations of Ai . Note that 
A a ^i and A a ^ do not commute unless a — 0, in which case the operators Tb,i an d 
To, 2 are equal to backward shifts and , respectively; in other words, the 
matrices 



c = [1 0] 



A 



0.1 



"0 


1 


0" 




"0 





1" 











, A . 2 = 
































(3.65) 



provide a commutative realization of Ai which is unique (up to unitary equivalence) 
by Theorem EHI Note also that the tuple (T a ,i,T a , 2 ) defined in (|3~H3) is never a 
contractive solution of the Gleason problem unless a = 0. 



Example 3.24. Consider the subspace 

4 Ai 



Ai = span 



4-AiA 2 '4-AiA 2 '4 
and define the operators T^i and T Q)2 on Ai by 

j3 + aaX 2 



Ai A 2 



cH(fe 



f 



4 - AiA 2 

where a s a fixed complex number and where 

4a + (3\i 



f 



7 + (1 — a)a\i 



AiA 2 



/(A) 



- 7^ 2 



4 - AiA 2 

is the generic element in Ai. Thus, /(0) = a and it is readily checked that 
aAiA 2 + 0\i + 7A 2 



/(A) - /(0) = 



4 - AiA 2 



= A 1 (T a4 /)(A) + A 2 (T a , 2 /)(A), 
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so the tuple (T 0l i,T 0) 2) solves the Gleason problem on Ai. As in the previous 
example, take the matrices 



C = [l 0] , 



where A a ^\ and A a ^ are the matrices of T aj i and T 0) 2 with respect to the basis 
i — - ^ ^— | of Ai and C is the matrix of the operator G : Ai — >• J 7 





"0 


l 

4 


0" 










1 

4 


A Q ,i = 













1 - a 










a 




















1 4— A1A2 ' 4— Ai A2 

defined in 



4— Ai A2 

For every choice of a, 



c(i - xUa,! - \ 2 A a<2 y 



Ai 



4-A1A2 



4-A1A2 



4-A1A2 



which realizes Ai as the range of the observability operator of a pair (C, A Q ). 
Different choices of a in (|3.fci4|l lead to non-equivalent realizations of Ai . Note that 
A a ^i and A a ^ never commute which is not surprising since Ai is not backward-shift 
invariant as has been established in Example 13. 191 

3.4. Applications of observability operators: the commutative setting. In 

this subsection we discuss applications of observability operators for the commuta- 
tive setting. This subsection parallels Subsection 12. 31 

For subspaces of TLy(kd) invariant under the forward shift operator-tuple Ma, 
we have the following analogue of the Beurling-Lax-Halmos-de Branges theorem 
due originally to Arveson [7] and McCullough- Trent (for the case of isometric 
inclusion); in fact, one can check that our proof, namely, the commutative adapta- 
tion of the proof of Theorem l2.14l follows that of |H] if one makes the substitution 
L = (0% , x*)* (where L is the key operator appearing in jHJ). In general, an op- 
erator between two Arveson spaces Tlu{kd) and Tty(kd) is said to be multiplier 
if intertwines the respective coordinate-function multipliers: 

0M Xj f = M Xj ef for all feH u {k d ). 

It is straightforward to see that a multiplier necessarily has the form 

Q = Mg: f(z)^d(z)-f(z) 

where 8(z) = J2nei^ @nZ n is a bounded, holomorphic C(U,y)-valued function on 

M d , but not all bounded, holomorphic, operator- valued functions on M d are mul- 
tipliers (see e.g. £Q). In case the multiplication operator has operator norm at 
most 1, we say that 9 is a contractive multiplier and belongs to the (commutative) 
multivariablc Schur-class Sd{U, y). Unlike the convention in the classical case, such 
a multiplier 9 is said to be inner if in addition Mg is a partial isometry. 

Theorem 3.25. (1) A Hilbert space M. is such that 

(a) M. is contractively contained in Jiyikd), 

(b) Ai is invariant under the Arveson-shift d-tuple M^, 

(c) the d-tuple 

M MX = {M MM M M ,X d ) where M M . X] := M Xj \m for j = 1, . . . , d 
is a row contraction: 

m m ,\A m mmT + ■ ■ • + m mm (m mm )* < i M , 

and 
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(d) (Mmx)* is strongly stable, i.e., 

W M M,x)* n ff M -^Oasn^^for all f e M 

nGZj : |n|=JV 

if and only if there is a coefficient Hilbert space IA and a contractive multi- 
plier 9 G Sd{U, y) so that M. = 8 ■ 7~iu(kd) with lifted norm 

\\0- f\\ M = \\Qf\\nu(k d ) 

where Q is the orthogonal projection onto (Ker Mg) C TLu{kd)- 
(2) The subspace A4 in part (1) above is isometrically contained in Tty(kd) if 
and only if the corresponding contractive multiplier 9 € SdiU,y) can be 
taken to be inner. 

Proof. The proof is a straightforward commutative adaptation of the proof of The- 
orem and hence will be left to the reader. We remark that, for the case where 
M is contained isometrically in ~Hy{kd), we are unable to obtain a representer 9 for 
which Mg is isometric but rather only a representer with Mg partially isometric. 
Indeed, one can check that the argument in the proof of Theorem 12. 141 breaks down 
because, for the case here, M Xj is only contractive rather than isometric. □ 

Remark 3.26. As observed in jS], from the function-theory point of view Theorem 
13.251 is not a true analogue of the classical Beurling-Lax theorem since the charac- 
terization of 9 is purely operator-theoretic with no information on the boundary 
behavior of the associated multiplier 9{z). This deficiency has now been remedied 
in the paper of Greene-Richter-Sundberg p?2] . 

The following is the analogue of Theorem 12.151 we omit the proof as it exactly 
parallels the proof of Theorem 12. 151 The result goes back to Drury |28| . 

Theorem 3.27. Suppose that T = (T l5 . . . , Td) is a commutative row-contractive 
operator-tuple with T* asymptotically stable and define the defect operator Z?x* o,nd 
the coefficient space y as in l|2.63[l . Then there is a subspace M. C 'Hyikd) invariant 
for the backward shift operator-tuple on Ti.y(kd) so that T is unitarily equivalent 
to PmM\\m- In particular, T has a Arveson-shift dilation unitarily equivalent to 
M\ on Hy{kd). 

As a corollary of this result one can arrive at the von Neumann inequality 

\\p{T u ...,T d )\\ < \\p(M Xl ,...,M Xd )\\ 

of Drury [2Hj and Arveson [Jj (see Remark 12.161 for the noncommutative case). 

Remark 3.28. The result in Theorem 13.271 is tied to the unit ball with associ- 
ated multivariable resolvent operator (/ — XiT* — • • • — AdTJ) -1 , associated defect 
operator Z?t* = (I — T\T* — ■ ■ ■ — X^T^) 1 / 2 , associated observability operator of 
the form Ofj T , = Z?t* (I — AiTj* — • • • — A^TJ) -1 and associated ambient kernel 
function fc(A, £) = 1/(1 — Aid — • • • — XdCd)- We mention that there has been a lot 
of work centering around other types of kernels and giving a model theory for other 
classes of operator-tuples by using appropriately modified observability-like opera- 
tors. Specifically, Miiller-Vasilescu [HHi for the commutative ball case with fc(A, £) = 
1/(1 — AiCi — ■ • • — \d(d) n \ Curto-Vasilescu [531 EU for the commutative polydisk 
case with fc(A,£) = (1/(1 — Aid) ■ ■ • (1 — ^dCd)) m , and Pott [Hj and Bhattacharyya- 
Sarkar JHj for the commutative case with fc(A, £) = 1/(1 — P(AiCi, • • • , ^d(d)) with 
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P equal to a "positively regular polynomial" . The most general form of results 
along this line is due to Ambrozie-Englis-Muller 0] and Arazy-Englis 0: given 
a positive-definite kernel k(X, C) on a domain D C C d and a d-tuple of operators 
T = (Ti , • • . , T c i) with Taylor spectrum contained in T> for which one can make 
sense of the defect operator D-r* := j:(T, T) and of the observability operator 

Od t ,,t*: x^D T ,k(X,T) 

(for example, if k(X, £) has no zeros inVxT) and T has Taylor spectrum contained 
in T>), then, under the assumption that Dt* > and that an additional stability 
condition on T* holds, 0.d t ».t* implements a unitary equivalence between T and 
PmM\\ m , where 

M = Ran Od t , , t » C U{k) ® y with y := Rlm"D T * , 

where M\ — (M\ 1 , . . . , M\ d ) is the operator-tuple of multiplication by the coordi- 
nate functions on TL(k)®y, and where Ai is invariant under each of , . . . , . 
The noncommutative case is not as well developed at this writing, but there is the 
paper of Popescu 0Hj which handles the case of a Cartesian product of noncommu- 
tative balls (and therefore including a noncommutative polydisk). We expect that 
many of the ideas of the present paper, including the interplay between the non- 
commutative and commutative settings and the connections with system theory, 
have some parallels in these other situations. 
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